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ABSTRACT 


The  KI^L^  Auger,  Is  and  2p  photoelectron  energies 

were  measured  for  a  series  of  phosphorus  and  sulphur  com¬ 
pounds  and  a  comparison  of  the  Auger  and  photoelectron  shifts 
was  made.  The  potential  model,  utilising  charges  from  both 
CNDO/2  and  extended  Hiickel  calculations,  was  used  to  attempt 
to  explain  the  shifts. 

On  comparing  the  Is  and  2p  photoelectron  shifts  a  good 
linear  correlation  was  obtained  with  the  Is  electron  showing 
the  larger  shift.  This  reflects  the  more  core-like  nature 
of  the  Is  electron  as  compared  to  the  2p  electron.  No  such 
correlation  was  obtained  on  comparison  of  the  Auger  electron 
shifts  with  the  photoelectron  shifts.  The  Auger  parameter, 
when  defined  properly,  related  the  Auger  electron  shift  to 
the  photoelectron  shift  in  a  convenient  way.  The  variations 
in  shift  could  be  rationalised  by  consideration  of  the  re¬ 
laxation  accompanying  the  loss  of  the  two  electrons  in  the 
final  two  hole  state.  It  was  found  that  relaxation  depended 
upon  the  ability  of  a  group  attached  to  the  central  atom 
to  release  electrons.  Many  factors  can  contribute  to  this 
and  reversals  of  the  expected  relaxation  trends  based  upon 
electronegativity  arguments  were  noted. 

The  group  shift  concept  was  also  applied  to  the  ex¬ 
perimental  data.  Inherent  in  the  group  shift  concept  is  the 
assumption  that  shifts  incurred  by  each  group  are  independent 
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of  each  other.  It  was  found  that  this  assumption  was  only 
valid  within  a  very  limited  series  of  compounds  at  best. 

The  ground,  relaxation  and  transition  potential  models 
were  tested  with  the  aid  of  two  semi-empirical  methods. 

The  CNDO/2  calculation  gave  reasonable  correlation  when 
using  fitted  parameters  though  the  parameters  obtained  did 
not  match  the  expected  theoretical  trends  too  well.  The 
EWMO  (Hiickel)  calculation  was  found  to  partition  charge 
differently  depending  upon  the  number  of  lone  pairs  on  the 
central  atom  and  hence  the  general  applicability  of  this 
method  was  severely  limited. 
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Chapter  1.  Introduction 


A.  General  Introduction 

In  1964,  Hagstrom,  Nordling  and  Siegbahn^  demonstrated 

that  high  resolution  3 -ray  spectrometers  could  be  used  to 

determine  the  binding  energies  of  inner  shell  electrons  in 

elements  of  low  atomic  number  and  that  the  binding  energy 

of  an  inner  shell  electron  from  a  particular  element  in 

different  chemical  environments  could  vary  significantly. 

Since  then  the  field  of  photoelectron  spectroscopy  has  been 

2 

shown  to  have  many  applications.  Photoelectron  spectroscopy 
is  in  essence  the  application  of  the  Einstein  Photoelectric 
Law: 

Er  =  hv  -  I  (1.1) 

which  relates  the  energy  of  the  incident  photon  (hv)  to 

the  ionisation  potential  (I)  and  the  kinetic  energy  of  the 

ejected  photoelectron  (E.J  ,  to  chemical  systems. 

l\ 

Following  initial  photoionisation  the  system  can  under¬ 
go  secondary  processes  to  achieve  a  lower  energy  state.  An 
electron  from  a  higher  energy  level  can  fill  the  vacancy 
releasing  the  excess  energy  either  as: 

(a)  a  photon  (X-ray  fluorescence)  or 

(b)  an  electron  (the  Auger  process) . 

3 

The  Auger  process  is  an  internal  rearrangement  of 
electrons  with  the  ejection  of  one  electron  carrying  the 
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excess  energy.  The  energy  of  the  ejected  electron  depends 
upon  the  energies  of  the  levels  involved.  As  is  the  case 
with  primary  photoionisation  the  Auger  electron  energy  can 
also  vary  significantly  with  the  chemical  environment. 
Understanding  how  these  shifts  arise  and  how  they  relate 
to  each  other  is  of  fundamental  importance. 

B.  Description  of  the  Various  Processes 

The  various  types  of  processes  which  can  occur  are 
illustrated  in  figure  1.1.  Note  that  both  X-ray  (K,  L,  M, . . ) 
and  orbital  (Is,  2s,  2p,  3s,  3p,._.)  notation  for  the  electron 
energy  levels  are  used  in  the  figure  to  demonstrate  the 
relationship. 

(i)  Photoionisation 

When  the  initial  state  of  a  molecule  or  atom  is  bom¬ 
barded  with  photons  of  a  characteristic  energy,  hv ,  a 
transition  occurs  in  which  the  final  state  is  an  ion  plus 
a  free  electron.  Conservation  of  energy  requires  that 

hv  =  ET,  +  E. 

K  1 

where  E,.  is  the  kinetic  energy  of  the  free  electron  and  E. 

K  l 

is  the  energy  of  formation  of  the  ion.  This  is  the  ionisation 

/ 

energy  for  removal  of  an  electron  from  the  molecule.  The 
ion  is  several  thousand  times  heavier  than  the  electron  and 
conservation  of  momentum  dictates  that  for  all  practical 
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Figure  1.1  Processes  involved  in  electron  spectroscopy. 

A.  Initial  state,  e.g.  Ne. 

B.  Photo  ionisation. 

C.  "Shake-up"  -  "Shake-off"  produces  satellite  structures. 

D.  Core  vacancy  produced  by  any  means  allows  E  and  F. 

E.  X-ray  fluorescence,  core  vacancy  filled  by  an  electron 
from  a  higher  level  and  a  photon  of  energy  released. 

F.  Auger  process,  core  vacancy  filled  by  an  electron  from  a 
higher  level  and  a  secondary  electron  is  released. 


■ 
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purposes  the  free  electron  takes  up  all  the  kinetic  energy. 

The  ionisation  energy  can  be  equated  to  the  difference 
in  the  total  energy  of  the  ion  state  and  that  of  the  ground 
state  of  the  molecule  or  atom: 

Ej8  =  EfS(N-l)  -  Eo(N)  (1.2) 

s 

where  =  ionisation  energy  of  electron  s 

Eq(N)  =  total  energy  of  the  ground  state  (the  initial 

state  of  the  system) 
s 

E^  (N-l)  =  total  energy  of  the  ion  which  is  formed  when 
electron  s  is  removed. 

Rigorous  calculation  of  the  appropriate  total  energies 

yields  the  ionisation  energy. 

Quantum  mechanically  the  probability  of  a  transition 

from  the  initial  ground  state  (Y")  to  the  final  state  (¥'  - 

ion  +  free  electron)  is  given  by  the  square  of  the  trans- 

4 

ition  moment  integral 

m  =  <  r  |  ep  |  r  >  d.3) 

where  p  is  the  dipole  moment  operator.  Application  of  the 
Born-Oppenheimer  approximation  separates  the  wave  functions 
into  a  product  of  electronic  and  nuclear  functions. 

Equation  1.3  then  becomes 

M  =  y?*"  (R)  V?v  1  (R)  dR.  JV*e"  (r;R)  |  Z  Pel^e*  (r;R)dr  (1.4) 
The  nuclear  function  has  been  further  split  into  rotational 


ri  • 
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and  vibrational  parts  and  since,  in  the  majority  of  cases 
rotational  structure  cannot  be  resolved,  the  former  is 
ignored . 

For  a  photoelectron  transition  to  be  allowed  the 
integrals  in  equation  1.4  must  be  non-zero.  The  final 
state  includes  a  free  electron  as  well  as  the  ion  and  hence 
there  is  always  a  non-zero  value  for  equation  1.4.  As  a 
consequence  all  one-electron  transitions  are  allowed. 

In  the  past  photoelectron  spectroscopy  has  been  divided 
(somewhat  artificially)  into  two  sections  depending  on  the 
photon  source  being  used.  One  major  branch  is  Photoelectron 
Spectroscopy  (PES)  or  Ultraviolet  Photoelectron  Spectroscopy 
(UPS) .  In  this  case  ultraviolet  radiation  provides  the 
photon  source  and  hence  this  spectroscopy  only  deals  with 
the  valence  region.  The  other  major  branch  is  Electron 
Spectroscopy  for  Chemical  Analysis  (ESCA)  or  X-ray  Photo¬ 
electron  Spectroscopy  (XPS) .  The  photon  source  in  this  case 
is  an  X-ray.  Depending  on  the  X-ray  being  used  and  the 
molecule  being  studied,  the  deeper  core  levels  of  the  atom 
can  be  probed  as  well  as  the  valence  shell  levels.  In 
this  thesis  only  the  techniques  appropriate  to  XPS  will 

5 

be  discussed. 

Because  the  core  electrons  do  not  take  part  in  bonding, 
it  was  originally  thought  that  they  would  be  of  little 
practical  interest  to  the  chemist.  It  is  now  known  that 
the  core  electron  binding  energy  values  are  affected  by  the 


' 
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valence  electron  distribution  and  yield  chemical  information 
upon  core  level  ionisation.  The  interaction  of  the  core 
electron  with  the  valence  shell  gives  rise  to  a  potential 
energy  contribution  to  the  total  ground  state  energy  of  the 
molecule.  Electron  withdrawing  substituents  would  reduce 
this  contribution  whereas  electron  donating  groups  would 
increase  the  contribution  and  hence  the  chemical  shift  from 
a  suitable  reference  gives  an  indication  of  the  chemical 
environment . 

(ii)  "Shake-up  and  Shake-off" 

In  addition  to  the  major  peak,  the  photoelectron  spectra 
of  the  core-levels  can  also  exhibit  many  satellite  struc¬ 
tures.  Energy-loss  peaks  occur  on  the  low  kinetic  energy 
side  of  the  main  peak.  They  arise  when  the  ejected  electron 
collides  with  molecules  in  the  chamber  thereby  losing  some 
kinetic  energy.  The  latter  contributions  are  pressure 
dependent  and  are  not  part  of  the  photoionisation  process. 

Upon  photoionisation  an  excited  ion  state  may  be 
formed  in  which  the  photoelectron  has  been  emitted  along 
with  the  simultaneous  excitation  of  an  outer  electron 
either  to  an  excited  bound  state  or  to  the  continuum.  The 
former  is  termed  "shake-up"  while  the  latter  is  termed 
"shake-off".  Upon  ionisation  the  resultant  ion  can  be  in 
one  of  a  number  of  states  with  the  ionisation  energy  of  the 
photoelectron  being  given  by  a  modified  form  of  equation  (1.2) 


I 
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Ej'k  =  E®'k(N-l)  -  Eo(N)  (1.5) 

When  k  =  0  the  ion  is  in  its  ground  state  and  this  gives 
rise  to  the  major  peak.  The  satellite  lines  are  defined  by 
k  =  1,2,  etc.  and  since  they  denote  excited  ion  states  the 
satellites  will  appear  on  the  low  kinetic  energy  side  of 
the  major  peak. 

Photoionisation  follows  dipole  selection  rules  (see 
equation  1.3),  that  is  AL  =  1  and  AS  =  0,  thus  the  allowed 
excited  ion  states  must  have  the  same  symmetry  as  the  ion 
in  its  ground  state. 

Satellites  also  appear  in  the  spectrum  due  to  the 
fact  that  the  photon  source  is  not  monochromatic.  The 
principal  satellites  are  observed  on  the  high  kinetic 
energy  side  of  the  main  peaks.  By  knowing  the  "diagram 
lines"  of  the  characteristic  X-ray  these  satellites  can 
be  easily  identified  and  removed  mathematically. 

(iii)  X-ray  Fluorescence  and  the  Auger  Process 

Following  creation  of  the  core-hole  state  of  the  ion, 
secondary  processes  can  occur  to  allow  the  system  to  achieve 
a  lower  energy  state.  The  two  major  decay  modes  are  X-ray 
fluorescence  and  the  Auger  process.  In  the  case  of  X-ray 
fluorescence  the  excess  energy  is  released  in  the  form  of 
a  photon  whereas  in  the  Auger  process  the  excess  energy 
is  given  to  an  emitted  electron. 
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For  the  lighter  elements  the  Auger  process  dominates.^ 
The  respective  yields  for  K-shell  fluorescence  and  Auger 
yields  are  given  by 


Pf+PA 


and  aK  =  P~i¥. 


(1.6) 


w  =  K-shell  fluorescence  yield 

Ja 

a.v  =  Auger  yield 

x\ 

=  transition  probability  for  X-ray  Fluorescence 
PA  =  transition  probability  for  the  Auger  process. 

For  P(Z  =  15)  w  =  0.061  and  a  =•  0.939.  Since  X-ray 

x\  X\ 

fluorescence  is  not  the  most  important  decay  mode  for  light 
elements  and  is  not  directly  measurable  with  our  equipment 
it  will  only  be  discussed  briefly  to  illustrate  the  dif¬ 
ferences  between  it  and  the  Auger  process. 

The  selection  rules  governing  X-ray  transitions  are 
given  by 


AL  =  ±1,  Aj  =  0,  ±1  with  j  =  0  j  =  0  (1.7) 


Thus  in  transitions  from  the  L  shell  to  a  K-hole  the  L2^Ka2^ 

and  L-,  (K  n  )  transitions  are  allowed  whereas  that  from 
3  al 

shell  is  forbidden. 

The  X-ray  line  emitted  corresponds  to  the  energy  dif¬ 
ference  between  the  two  levels  involved  in  the  X-ray 
transition,  and  to  the  first  approximation  the  shifts 
for  the  core  levels  of  a  molecule  are  very  similar. 


Thus 


. 
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shifts  in  the  X-ray  line  in  a  series  of  molecules  should 
be  very  similar.  This  limits  the  use  of  X-ray  emission  for 
the  study  of  the  chemical  environment  however,  as  will  be 
seen  later,  there  is  a  measurable  difference  in  core 
electron  chemical  shifts  between  the  Is  and  2p  levels 

p 

(1.3  eV  for  PF^  compared  to  PH^) . 

The  energy  of  the  Auger  electron  is  given  by  the 
difference  between  the  total  energy  of  the  initial  hole 
state  and  that  of  the  two  hole  state: 


eau(xyz)  em+(x)  em++(yz) 


(1.8) 


Ean (XYZ )  is  the  kinetic  energy  of  the  emitted  XYZ  Auger 
electron,  EM+(X)  the  total  energy  of  the  initial  hole 
state,  a  singly  ionised  species  with  a  hole  in  level  X, 
EM++(YZ)  the  total  energy  of  the  final  two-hole  state, 
a  doubly  ionised  species  with  holes  in  levels  Y  and  Z. 
Thus  in  figure  IF  the  emitted  Auger  electron  would  be 
designated  KL^L^ .  The  nomenclature  follows  the  j-j 
coupling  scheme  which  will  be  described  later. 

The  Auger  energy  can  also  be  expressed  in  terms  of 

9 

the  binding  energies  of  the  levels  involved: 


eau(xyz)  =  eb(x)  ~  VY)  "  EB(Z)  “  R 


(1.9) 


R  includes  terms  which  further  reduce  the  Auger  energy  due 
to  differences  in  orbital  energy  values  once  an  electron 
is  removed.  This  form  allows  the  Auger  energy  to  be  put  in 


’ 


terms  of  measurable  quantities  and  provides  a  basis  upon 
which  semi-empirical  methods  of  calculation  can  be  applied. 
The  indistinguishability  of  electrons  does  not  permit 
distinction  of  whether  electron  Y  or  Z  is  emitted  as  the 
Auger  electron. 

Unlike  X-ray  fluorescence,  the  Auger  process  can 
involve  the  shell.  The  selection  rules  are:^ 


AL  =  AS  =  AJ  =  0  ,  parity  unchanged 
and  the  transition  probability  is  given  by: 


-  TT 


(r2)lfrjr^l*i(rl)  Vr2>  (1* 


The  Auger  transition  arises  from  a  coulombic  rearrangement 
due  to  the  interaction  of  the  two  electrons  involved. ^ 

<Jk  and  tJk  are  the  single  electron  wavefunctions  describing 
the  initial  hole  state  and  and  are  the  wavefunctions 
describing  the  final  state,  which  includes  the  continuum 
electron.  The  wavefunctions  must  be  correctly  anti-sym¬ 
metrised.  The  continuum  electron  can  be  treated  as  a  plane 
wave  and  hence  its  orbital  angular  momentum  is  unrestricted 
thus  allowing  many  transitions. 

The  second  row  (Si-Ar)  KLL  Auger  spectra  provides  a 
suitable  system  for  the  study  of  Auger  electron  chemical 
shifts.  Figure  1.2  shows  the  Auger  KLL  spectra  of  Ar,  PH^ 
and  PF^ .  The  spectra  are  atomic  in  nature  and  show  very 
similar  features.  Auger  spectra  which  involve  the  valence- 
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The 


Figure  1.2  KLL  Auger  spectra  of  Ar,  PH^  and  PF^. 

initial  hole  state  was  created  with  the  Ti  K  X-ray 

a 

(4510.84  eV) .  The  major  peaks  on  the  Argon  spectrum  are 
labelled  in  both  the  j-j  and  L-S  coupling  schemes.  The 
'  1  *  denotes  satellites  thought  to  be  associated  with  M 
electron  "shake-off". 
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shell,  the  KLL  Auger  spectra  of  the  first  row  for  example,"^ 
are  more  complex  because  the  L  shell  atomic  orbitals  in¬ 
volved  are  now  shared  with  other  atoms.  Hence  KLL  Auger 
spectra  of  these  elements  are  unsuitable  for  the  study  of 
chemical  shift,  however,  complimentary  information  to  UPS 
can  be  derived  therefrom. 

The  number  of  lines  in  the  low  Z  region  of  the  periodic 

table  can  be  described  by  L-S  coupling.  In  this  case  the 

spin  angular  momenta,  s>  of  each  of  the  electrons  involved 

couple  to  give  the  total  spin  orbital  momentum,  S,  of  the 

system  and  similarly  the  orbital  angular  momenta,  £,  of 

each  electron  couple  separately  to  give  the  total  orbital 

angular  momentum,  L,  of  the  system.  The  total  spin  and 

orbital  momenta  couple  to  give  the  total  angular  momentum, 

J,  of  the  system.  The  resultant  state  is  designated  as 
2s+l 

L  where  2S+1  is  called  the  multiplicity,  J  is  the 
J 

total  angular  momentum  and  L,  the  angular  momentum  is 
designated  by  the  letters  S,  P,  D,  F,  etc.,  for  L=0,  1,  2,  3 
respectively.  Six  final  states  are  possible  and  these  are 
shown  in  table  1.1. 

The  KLL  Auger  spectrum  shows  only  five  major  lines 

3  2  4 

in  the  low  Z  region  since  the  P  state  from  the  2s  2p 

configuration  is  forbidden  because  parity  is  not  conserved. 

The  initial  K-hole  state  (ls1/2^  in  the  KLL  Au9er  process 

has  even  parity.  The  final  state  consists  of  a  doubly- 

2  4  3 

ionised  core  hole,  which  in  this  case  is  the  2s  2p  (  P) 


, 


' 
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TABLE  1.1 


CONFIGURATION  STATES 

2s°2p6  1S() 

2sl2P5  3p0 ,1,2'  S 

2s22P4  S''  3p0,l,2'‘S 


configuration  having  even  parity  and  an  emitted  electron. 

The  electron,  a  continuum  wave,  must  be  in  the  P  state 

(i.e.  of  angular  momentum  L  =  1)  in  order  to  conserve  the 

total  angular  momentum  of  the  system.  Hence  the  final 

state  of  the  total  system  will  have  odd  parity  thereby 
3 

making  the  P  state  forbidden. 

In  the  high  Z  region  j-j  coupling  dominates.  In 
this  scheme  the  s  and  Z  momenta  of  each  electron  couple 
to  form  a  total  angular  momentum,  j,  for  each  electron. 

Thus  each  electron  is  described  by  a  definite  j  value  which 
must  be  specified  before  electron-electron  coupling  is 
invoked.  This  gives  the  23.jy2  (L^)  ,  2pjy2(L2)  an^  ^P3/2^L3^ 
electrons,  and  along  with  the  initial  hole  state,  ^S^y2 (K) / 
gives  rise  to  six  possible  Auger  lines,  namely;  KL^L^ , 

KL^L2  ,  YJj^L2  1  kl-lL3  9  KI^L^  and  KXj^Li^  . 


By  following  the  total  J  values,  the  j-j  coupling 

scheme  can  be  correlated  with  the  L-S  scheme.  This  is 

illustrated  in  figure  1.3. 10  In  the  intermediate  coupling 

3 

scheme  nine  Auger  lines  are  now  possible.  The  PQ  and 
3  2  4 

P2  lines  of  the  2s  2p  configuration  are  able  to  mix  with 
the  1S0  and  1D2  lines  of  the  same  configuration  res¬ 
pectively  via  the  same  J  quantum  number  and  parity  is  no 

3  2  4 

longer  violated.  The  P1  state  of  the  2s  2p  configuration 

2  4 

is  still  forbidden  since  there  is  no  state  of  2s  2p 
configuration  with  the  same  J  quantum  number  for  it  to  mix 


with . 


. 


Relative  energies 
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Intermediate 


Figure  1.3  KLL  line  positions  as  a  function  of  atomic 
number.  (Reproduced  with  permission  from  Reference  10.) 
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The  major  lines  are  marked  on  the  Ar  spectrum  (Figure 
1.2)  using  both  coupling  schemes.  In  this  region  of  Z, 
the  L-S  coupling  still  dominates  and  only  the  five  major 
lines  are  observed.  Other  peaks  are  also  present,  the 
majority  of  these  peaks  arise  from  the  emission  of  an  Auger 
electron  associated  with  a  more  highly  excited  initial  K 
hole  state.  Thus  Auger  electrons  emitted  by  such  excited 
states  will  have  a  lower  kinetic  energy. 

The  excited  initial  K-hole  state  can  be  produced  by 
"shake-up"  or  "shake-off"  accompanying  the  initial  photo¬ 
ionisation.  For  the  second-row  elements  the  probability 
of  outer  shell  excitation  accompanying  K  ionisation  is 

approximately  4%  from  the  L  shell  and  20%  from  the  M 
12 

shell.  Peaks  associated  with  concomitant  M  shell  ex¬ 
citation  should  appear  on  the  low  kinetic  energy  side  of 

each  peak  and  should  possess  similar  relative  energies  and 

13  14 

intensities  as  the  major  Auger  lines,  '  since  the  M 
shell  electrons  are  not  involved  in  the  KLL  Auger  process 
and  should  not  interact  strongly  with  the  L  shell  electrons. 
The  lines  marked  '1'  in  Figure  L  2  are  thought  to  be  satellite 
lines  which  are  associated  with  M  electron  "shake-off"  and 
can  therefore  be  designated  as  KM-LLM  lines.  It  can  be 
seen  that  they  are  equidistant  from  their  respective  ground 
state  Auger  peaks.  Comparing  the  spectra  of  PH^  and  PF^ 
reveals  that  the  energies  of  these  features  are  chemically 

TO  A 

dependent.  '  Peaks  arising  from  states  in  which  the  L 
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shell  is  excited  are  less  intense  than  those  in  which  the 

M  shell  is  excited.  These  satellites  should  be  irregular 

both  in  position  and  intensity  since  the  excited  electron 

arises  from  the  same  shell  as  that  partaking  in  the  Auger 

14 

process,  and  strong  interactions  would  be  expected.  Ex¬ 
citation  of  a  bound  electron  upon  ejection  of  the  Auger 
electron  can  also  contribute  to  the  satellite  structure. 

This  is  known  as  the  double-Auger  process  and  should  be 
difficult  to  see  because  the  energy  is  divided  between  two 
electrons . ^ 

Excitation  of  the  Auger  spectrum  can  be  effected  by 
different  means;  the  only  essential  feature  is  the  creation 
of  an  ion  with  an  inner  shell  hole.  X-ray  excitation  used 
herein  generally  yields  a  simpler  spectrum  because  production 
of  highly  excited  ions  is  limited.  Electron  beam  excitation 
frequently  produces  highly  excited  hole-state  ions,  hence 
satellite  features  are  frequently  more  prominent. 

C.  Models  for  Chemical  Shift 

It  is  not  feasible  or  practical  to  calculate  initial 
and  final  state  energies  for  all  molecules  studied  in  order 
to  describe  the  chemical  shift  in  both  photoelectron  or 
Auger  spectroscopy.  Simple  models  can  greatly  extend  the 
utility  of  the  data  obtained  and  we  discuss  here  the  basis 
for  the  models  and  the  concepts  involved. 
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(i)  Koopman's  Theorem  versus  AE 


.  The  Concept  of  Relaxation 


SCF 

The  binding  (or  ionisation)  energy  is  properly  defined 
by  the  difference  in  the  total  energy  of  the  ground  state 
and  that  of  the  final  state  (equation  1.2).  The  total 
energy  of  the  ground  state  can  be  obtained  by  use  of  a  full 
self  consistent  field  calculation  within  the  Hartree-Fock  scheme, 
as  can  that  of  the  final  state  with  the  appropriate  core 
hole.  The  difference  in  the  total  energies  will  give  the 
binding  energy  of  the  electron.  This  is  known  as  the  AEgCF 
method.^  The  method  can  prove  expensive  and  within  itself 
be  subject  to  various  levels  of  approximation  and  sophis¬ 
tication  such  as  the  extent  and  type  of  basis  set  and  the 
inclusion  of  configuration  interaction. 

The  most  drastic  simplification  that  can  be  made  to  the 
scheme  is  to  use  the  same  wave  functions  for  the  passive 
orbitals  in  the  ionised  state  as  in  the  ground-state.  This 
is  the  "frozen  orbital  approximation". 

The  Hamiltonian,  in  the  Hartree-Fock  scheme,  for  an 
N-electron  system  with  spin-orbitals  Y.  which  are  solutions 
of  the  Hartree-Fock  equation  H¥ .  =  is  given  by 

N 

H  =  h  +  Z  (J.-K.)  (1.11) 

j=l  3  3 

12  Za 

where  h  is  the  one-electron  operator  (-■ -  Z  - — )  and  J. 

2  a  r la  3 
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Kj  are  the  coulomb  and  exchange  operators. 

The  total  energy  for  the  N-electron  system  is 


N  N 

=  £  h.  +  ~  Z  Z 

i=l  1  2  i=l  j=l 


(J. .-K. . ) 
1J  13 


(1.12) 


where  =  <4'±  (1)  V  ^  (2)  |  l/r12  |  V  ±  (1)  f.(2)> 

Kij  =  (2)  |l/r12|l'j  (1)  Ti(2)> 

and  for  the  N-l  system,  with  electron  £  missing 


N  N 

E  (N-l ,  £ )  =  E  h.  +  i  E  E  (J.  .-K.  .)  (1.13) 

iftl  1  13  13 


The  difference  is  simply 


E  (N) 


E  (N-l ,  £  ) 


i  N 

h*  +  2  E 


(1.14) 


but  this  is  simply  the  eigenvalue  for  electron  i(e  )  in 

/v* 

the  original  Hartree-Fock  equation 


E(N)  -  E  (N-l , £  )  =  (1.15) 

15  17 

This  result  is  known  as  Koopman 1 s  theorem  '  which  states 
that  the  binding  energy  is  taken  to  be  the  negative  of  the 
eigenvalue  of  the  Hartree-Fock  equation.  Only  one  calculation 
is  required. 


' 
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The  major  flaw  in  the  Koopman ' s  Theorem  approach  is 
that  it  is  a  frozen  orbital  approach  and  takes  no  account 
of  electronic  reorganisation  which  accompanies  photoion¬ 
isation.  This  would  not  be  a  problem  if  reorganisation 
energies  were  constant  for  each  element,  however  it  is 
apparent  that  the  contributions  vary  significantly  with 
structure  and  valence.  A  method  of  correcting  the  Koopman' s 
Theorem  value  to  the  full  ab  initio  AE_pt;,  standard  would 
provide  a  cost  and  time  advantage,  avoid  problems  that  can 
possibly  occur  with  hole-state  calculations^  as  well  as 
providing  a  deeper  conceptual  understanding  of  the  nature  of 

the  processes  involved.  Such  an  approach  was  described  by 

1 8 

Hedin  and  Johannson. 

The  Hamiltonian  for  the  ground  state  was  written  as: 


H  =  h  +  I  V.  =  h  +  V 
k  k 


(1.16) 


h  is  the  one-electron  operator,  =  J^-K^,  the  two-electron 
coulomb  and  exchange  operators  from  orbital  k.  For  the 
hole  state  the  Hamiltonian  was  written  as: 


/\  *  * 
H  =  h  +  V 


(1.17) 


where  V  =  V-V,  +  V  ,  the  *  indicating  a  hole  state.  The 

k  p 

potential  describing  the  hole  state  is  the  original  potential 
minus  that  due  to  the  removed  electron  from  orbital  k  plus 
a  polarisation  potential  due  to  the  polarisation 


■ 
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of  the  orbitals  upon  removal  of  electron  £. 


V 

P 


I 

k?l 


<w 


(1.18) 


This  is  the  sum  of  the  differences  in  the  potential  of  the 

passive  orbitals  before  and  after  ionisation. 

18 

Hedin  and  Johannson  showed  that 


£a  +  I 


<l\v  \l> 

1  p1 


E  <6 i | H* 


-H-6e  .  I  <5i> 

l 1 


(1.19) 


the  6  signifying  the  difference  between  the  initial  and 
final  values.  Expression  (1.19)  is  exact,  no  approximations 
have  been  made.  The  final  term  on  the  right-hand  is  small 
compared  to  the  polarisation  (or  relaxation)  term  and  so 
the  relaxed  orbital  energy  becomes 

E*  =  ex  +  j  <i|Vp|t>  (1.20) 

This  result  provides  a  good  basis  for  further  studies, 
especially  in  the  formulation  of  semi-empirical  potential 
models . 

(ii)  Correlation  of  Potential  with  Binding-Energy  Shifts 
19  20 

Basch  and  Schwartz  showed,  almost  simultaneously, 
that  the  change  in  orbital  energy  paralleled  the  change  in 
potential  at  the  nucleus  on  which  the  core  orbital  resided. 

The  potential  at  nucleus  A  arising  from  doubly -occupied 


23 


orbitals  k  and  other  nuclei,  Z  ,  is  given  by 

j3 


=  -21  <k|l/r1A|k>  +  £  Zd/R 


B^A 


B'  AB 


(1.21) 


In  the  Hartree-Fock-SCF  theory  the  orbital  energy  is  given 
by 

£k  =  <kl’|Vllk>  +  <kl^  "  ZB/rlBlk>  + 


Z 

i^k 


K,  J 
k£' 


+  J 


kk 


(1.22) 


This  is  the  same  as  equation  (1.14)  except  that  the  sum¬ 
mations  are  performed  over  doubly-occupied  orbitals. 

Basch  collapsed  orbital  k  onto  nucleus  A  for  all 
terms  except  those  giving  zero  or  indeterminate  results. 

and  Z  K,  0  were  kept. 


Thus  <k|-|v^|k>,  Jkk 


<k 


Z  ZB/rlB 

D 


k> 


<k|-VrlAlk>  +  <k,B^A"VrlBlk> 

<kl"ZA/rlA>k>  +  Z_ZB/RAB 


(r^B  -  R  ,  the  internuclear  distance  since  electron  1  is 
centered  on  A. ) 


Z 

Jt^k 


2  Z 
SL/k 


<k  (2)  l  (1)  |  l/r1A|k(2)  t  (1)  > 


2  £  <k  |  k><!  1 1/r,  ,  |  !■>  =  2  £  <i|l/r  la|A> 
Iji  k  1  A*k  1A 


:*r  '  *v|-  \*>WM 


» 
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Hence  equation  (1.22)  became 


£k  “  VA  +  <k!-^llk>  +  <k|-VrlAlk>  +  Jkk'^k  Kk*  (1 


where  V  =  2  Z  < «.  X/r  Jl>  +  z  -Z  /R 

A  J#k  iA  B^A  B  AB 

and  finally  e  ^  AV 

J»>  r\ 


on  the  assumption  that  the  change  in  the  other  terms  in 

equation  (1.23)  can  be  taken  to  be  negligible.  Calcula- 
19 

tions  have  indicated  that  this  assumption  is  valid. 

Comparison  of  with  $  (equation  1.21)  shows  it  to 

be  basically  an  "external  potential",  with  AV  being  the 

difference  in  potential  "felt"  by  the  atom. 

20 

Schwartz  analysed  the  system  in  a  slightly  dif¬ 
ferent  manner.  The  Is  orbital  was  separated  from  the 
others  which  had  previously  been  transformed  to  localised 
molecular  orbitals,  L.  The  interaction  of  the  Is  orbital 
with  orbitals  L  in  the  two  electron  integrals  was  further 
split  into  a  local  contribution,  (molecular  orbitals 
connected  to  atom  A)  and  a  distant  contribution,  (those 
not  connected  to  atom  A).  Equation  1.22  became 


eis  ■  <lsl-|vrVriAlls>  +  Jisis-B^A<lsl  VriBlls>  + 


S  (2Jlsi-KlSi>  +  ,Ia1,t<2Jlsj-Klsj) 
i=loc  j=dist 


23) 


(1.24) 


■  tt> 


25 


This  form  allows  the  following  approximations  to  be  made. 

The  Is  orbital  is  highly  localised  (becoming  essentially 
a  delta  function)  and  so  the  distant  exchange  integrals 
(K^sj)  should  be  approximately  zero.  Furthermore,  the 
distant  coulomb  integrals  should  have  the  same  order  as  the 
electron-nuclear  attraction  integral  of  the  distant  for 
nucleus  A  and  the  attraction  of  the  Is  density  to  other 
nuclei  can  be  approximated  by  a  point  charge: 


<1SlZB/rlBlls>  =  VRAB 


Finally,  the  localised  internal  energy  of  Is 


4T=  <isi-lvi  -  vriAiis>  +  jisis 


should  be  insensitive  to  environment,  leaving 


-e 


ext 


int 

£ls  cls 


=  -2  E 


<L  .  !  1/r ,  A  I  L  .  >  +  E  Z  /R 


j=dist  J'  '  1A  j  B^A  B'  AB 


.Z  (2Jlsi  Klsi) 
i=loc 


(1.25) 


comparing  this  with  $ext 


$  =  $  -  $ 
ext  int 


-  -2  E.  <Lj  |  1/riAl  Lj>  +  ZB/RAB 

j=dist  J 


B^A 


. 


■ 
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2  E  < 
i=loc 


Liii/riAiLi> 


(1.26) 


The  final  expressions  on  the  right-hand  side  of  equations 
(1.25)  and  (1.26)  are  not  identical  to  each  other  but  should 
be  similar  and  their  changes  with  the  environment  should 
be  very  nearly  the  same,  thus: 


A<-els>  = 


A(W 


and  by  application  of  Koopman * s  Theorem 


AEb  »  A 


<*ext> 


Furthermore,  the  core  orbitals  at  the  other  nuclei  are 
very  localised  and  screen  their  respective  nuclei  and  so 
only  the  valence  orbitals  need  be  considered: 


AEB  =  A($val> 


(1-27) 


where  $ 


val 


=  -2  Z  <£|l/r1A|£>  +  I  Zn/RaT3  with  Zn  being 


B'  AB 


B 


£=val  B^A 

the  effective  reduced  nuclear  charge.  This  result  allows 

semi-empirical  models  utilising  such  methods  as  the  CNDO 

21 

formalisation  to  be  used. 


The  potential  model  can  also  be  expressed  in  the 
following  manner  for  an  N-electron  system 


F  -  kN  aN  +  VN  +  £ 
EB  -  ki  qA  +  VA  +  *■ 


(1.28) 


this  expression,  which  will  be  termed  the  "atom-charge 


'  .  '  •  ^fjfj  I  1 1 


. 
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model",  was  originally  based  upon  purely  electrostatic 
considerations.  It  can  also  be  shown  to  follow  directly 
from  the  Hartree-Fock  orbital  energy  using  arguments 

1  Q  Tf] 

analogous  to  those  of  Schwartz  and  Basch.  ' 


q^  =  charge  on  atom  A 


Va  =  2  q/R  ;  q  =  charge  on  atom  B 

A  B^A  B  AB  B 


Rab  =  internuclear  distance  between  atoms 
A  and  B 


.N 


Thus  is  the  potential  on  atom  A  due  to  the  charges  at  all 
other  atoms. 


k^  =  <ij  1 1/r^j  |  ij> ;  j  =  valence  electron. 

This  is  the  coulomb  repulsion  integral  between  the  core  and 
the  valence  electron. 

The  value  i  is  simply  a  constant  depending  on  the 
reference  level. 

By  invoking  the  point-charge  approximation  the  atom 
charge  model  reduces  to  the  potential  model  of  Schwartz  or 
Basch.  In  this  case 


ki  =  <  j  1 1/rJ  j>  =  C/n 

The  nuclear  attraction  integral  of  the  valence  orbital  at 
A  for  nucleus  A,  which  in  turn  can  be  equated  to  the  Slater 
orbital  exponent  for  the  valence  shell  (C)  divided  by  the 
principal  quantum  number  (n)  of  that  shell  provided  Slater 


' 


. 

■  ■ 
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type  orbitals  are  used.  Both  the  calculated  integral  a 

24 

and  the  Slater  exponent  have  been  employed  in  the 
evaluation  of  this  model. 

The  value  for  k_  can  also  be  determined  empirically 
by  a  least-squares  fit  on  experimental  data  in  conjunction 
with  the  charges  and  potentials  obtained  from  the  semi- 
empirical  calculation.  A  separate  k^  can  thus  be  obtained 
for  each  core-level. 

The  above  models  have  assumed  the  validity  of  Koopman's 
Theorem.  This  is  satisfactory  for  predicting  shifts  in 
carefully  selected  series  of  compounds  where  the  relaxation 
is  similar.  This  model  is  known  as  the  Ground  Potential 
Model  (GPM) .  Further  improvement  would  allow  for  variation 
in  relaxation  and  is  known  as  the  Relaxation  Potential 
Model  (RPM).23 

The  relaxation  potential  model  follows  directly  from 
18 

equation  (1.20).  Re-expanding  this  equation: 

Ej,  =  <  2. 1  h+v|  «■>  +  |<Jl|  r. 

k^  & 

=  i<£|h+v|£>  +  ^Ih*  + 

-  +  (1-29) 

■k 

since  by  definition  <£[v^|£>  =  < |  l/r^  I  =  0 >  is  "the 

hole  state  orbital  energy.  The  result  confirms  an  empirical 

25 

rule  proposed  by  Liberman.  Equating  the  change  in  orbital 
energies  to  changes  in  potential  leads  to 


- 


I  ■ ■■ . 
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E*  ■  b  VA  +  V  =  -  eb  (1-30> 

where  is  the  shift  in  the  binding  energy  of  the  core¬ 

electron.  One  further  approximation  must  be  made  if  semi- 

•  *  T  * 

empirical  approaches  are  to  be  used.  To  calculate  V 

ri 

the  concept  of  equivalent  cores  must  be  invoked.  The  idea 

was  originally  applied  to  the  estimation  of  binding  energy 

2  6 

shifts  by  Jolly  using  thermochemical  data.  However  the 

27 

concept  had  been  well  established  for  analogous  systems. 

The  Is  electrons  effectively  shield  the  other  electrons 
from  one  unit  of  nuclear  charge.  In  the  (Is)  hole  state 
the  valence  electrons  belonging  to  a  nucleus  of  charge  Q 
are  subjected  to  (Q-l)  units  of  charge.  This  behaviour  can 
be  mimicked  by  replacing  the  parameters  of  the  ionised 
atom  by  those  of  the  atom  with  one  additional  unit  of 
atomic  charge,  ie  Z+l.  Increasing  the  net  molecular  charge 
by  one  unit  equates  the  number  of  electrons  in  the  two 
species  and  preserves  the  closed  shell  structure.  Thus 
the  M(Z)Xn  core  hole  ion  is  approximated  by  the  M(Z+l)Xn+ 
ion.  The  same  number  of  valence  electrons  in  the  ion  are 
subjected  to  a  nucleus  of  total  charge  (Q+l)  shielded  by 
two  Is  electrons  and  hence  are  under  the  influence  of  (Q-l) 
units  of  charge  as  in  the  case  of  the  ion  with  a  (Is)  hole. 
The  assumptions  made  in  the  equivalent-core  approximation 
are  on  the  same  level  as  those  leading  to  the  potential 
model^^C  (compare  equation  1.27)  and  so  equation  1.30 


2  ■ 


■ 


. 
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finally  becomes 


aeb  =  "I  (ava(z>  +  ava(z+1)) 


U. 31) 


(iii)  Models  Involving  Transition  Operator  Formalism 


The  RPM  requires  that  two  calculations  be  performed. 

2  8 

The  Transition  Potential  Model  (TPM)  has  been  suggested 
as  a  method  of  avoiding  the  need  for  two  calculations. 

The  TPM  follows  from  the  Transition  Operator  Method  (TOM) 


29 


and  before  discussing  the  specific  application  it  is  nec¬ 
essary  to  consider  some  of  the  formalism  behind  TOM. 

TOM  maintains  the  one-electron  picture  at  AESCF  ac" 
curacy  by  optimising  the  mean  energy  of  the  initial  and 
final  states  with  respect  to  a  set  of  common  spin  orbitals. 
The  shift  in  binding  energy  is  obtained  from  the  eigenvalue 
of  the  ionised  spin  orbital  (i)  using  the  transition  Fock 
operator  which  is  defined  as 


fT (1)  =  hx  + 


<jT(2)  |  |  jT(2)>  +  i-  <ii(2)|  |ii(2)>  (1.32) 


•  T 


.T 


(where  <a(2)||a(2)>  q(l)  =  / a  (2)  l/r^  ^”P12^  a^)  q(l) 
and  associated  with  the  following  eigenvalue  problem 

fT  (1)  qT  (1 )  =  eT  qT(l)  (1.33) 


The  coulomb-exchange  term  associated  with  the  i-th  spin 
orbital  is  multiplied  by  occupation  number  to  simulate 
the  mean  initial  and  final  state. 


31 


Second-order  perturbation  analysis  has  shown  that  the 

"transition"  eigenvalue  differs  from  the  Koopman's  Theorem 

►  2  9 

value  in  the  same  way  as  does  the  AE^^,  ionisation  energy: 

“AESCF^  ~  ei  ~  ei+  2  |  <  j  i  |  |  ai>l2/  (e  -e  • )  (1.34) 

j  ra  J 


where  a  is  a  virtual  orbital.  Further  analysis  has  led  to 

T 

a  variational  argument  to  justify  F  and  it  has  also  been 
shown  that  the  "transition"  eigenvalue  is  equal  to  the 

AESCF  energy  through  to  third-order  in  perturbation 

t  .  30 

analysis . 

With  the  justification  of  TOM  it  is  now  convenient  to 

extend  the  transition  concept  to  the  potential  model.  The 

ionised  atom  in  the  molecule  of  interest  is  replaced  by  a 

*  1 

pseudo-atom  of  effective  charge,  Z^  =  Z^  +  ^ •  The  para- 

V 

meters  used  for  this  atom  in  semi-empirical  calculations 

are  interpolated  between  those  appropriate  to  the  Z^  and 

Z_  +  1  values.  Put  in  the  form  of  the  atom  charge  model 
A 

(equation  1.28),  the  TPM  can  be  expressed  as 


AE  =  kN'N-X 
B 


N,N-1  , 
+ 


V  N'N"1  + 
A 


(1.35) 


The  superscript  (N,N-1)  denotes  (transition)  values  as¬ 
sociated  with  a  system  of  going  from  N  to  N-l  electrons. 
Comparison  of  the  TPM  with  the  GPM  should  yield  an  estimate 
of  the  relaxation  energy. 


' 


■ 
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(iv)  Extension  of  Potential  Models  to  Auger  Spectroscopy 


The  Auger  KLL'  electron  kinetic  energy  is  given  by  the 
difference  in  the  total  energy  of  the  initial  state,  a 
species  with  a  K-shell  hole,  and  the  total  energy  of  the 
final  state  which  contains  2  L  shell  holes. 


ea(kll' ) 


(N-2)  (N-2) 

K  "  ELL' 


(1.36) 


In  terms  of  the  binding  energies  of  the  levels  involved, 
the  Auger  process  can  be  expressed  as: 


Ea(KLL')  =  Efi(K)  -  eb<l)  "  eB(Ij,)  (1*37) 


The  bar  indicates  that  electron  L'  arises  from  a  hole  state. 

Equation  1.37  is  equivalent  to  equation  1.9. 

Both  Liberman  and  transition  operator  approaches  can 

31 

be  applied  to  these  equations.  Adams  applied  the  Liberman 
expression  to  the  Auger  process  for  some  chlorine  containing 
species.  Equation  1.36  was  modified  to  become  a  " two-step" 
process . 


ea(kll' ) 


(E 


N-l 

K 


(E 


N-l 

LL' 


=  EB(K)  -  Eb(LL‘) 


(1.38) 


E  (LL1)  can  be  called  a  "two  electron  binding  energy". 

B 

Applying  the  Liberman  expression  to  equation  (1.38)  gave: 
ea  =  -^-[e  (K) N  +  TdO1^”1]  +  i[e(L)N  +  e(L)N  2]  + 


P'L 


< 

"m'r 

’ 


' 
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|[e(L')N  +  F(L')N~2]  (1.39) 

and  to  equation  1.37  ("three-step"  process) 

Ea  =  -|[e(K)N  +  F(K)N_1]  +  |[e(L)N  +  e(L)N-1]  + 

i[e  (L'  )N_1  +  e  (L'  )N_2]  (1.40) 

Both  of  these  expressions  were  tested  against  cal- 

otr 

culations  and  it  was  found  that  the  three-step  process  gave 

31 

the  better  result.  Adams  applied  the  "potential  at  the 
nucleus"  approximation,  the  equivalent  core  approximation 
and  made  the  further  (dubious)  approximation  (which  is 
inherent  in  the  "potential  at  the  nucleus"  approach)  that 
core-level  shifts  were  the  same,  obtained  the  following  RPM 
expressions : 

AEA  "  I[AVN(Z)  ”  AVn(Z+1)  +  2AVN<Z+2)]  (1.41) 

for  the  two-step  process  and 

AEa  *  |[AV  (Z+l)  +  AVn(Z+2)]  (1.42) 

for  the  three-step  process.  Equation  1.42  is  also  the 
negative  of  the  RPM  expression  for  the  binding-energy  shift 
of  a  core-electron  leaving  an  ion  with  a  single-core  hole. 
Thus  by  utilizing  equation  1.37  the  Auger  energy  shift  can 
be  written  as 

AEa  =  AE^*  +  |(AV(Z+1)  +  AV  (Z+2)  ] 


(1.43) 


■ 


- 


. 
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where  AE 


Lr*K 


AEb(K)  -  A  Eg  CL) 


X-ray 

The  Is  level  and  2p  level  interact  differently  with 
the  valence  shell  due  to  their  different  proximities.  This 
manifests  itself  in  a  larger  shift  for  the  Is  electron  than 

g 

for  the  2p.  The  "potential  at  the  nucleus"  model  should 
reflect  the  Is  shift  better  than  the  2p,  however  the 
potential  part  of  equation  1.43  refers  to  an  electron  from 
the  2p  shell.  A  measure  of  the  interaction  of  the  core¬ 
levels  with  the  valence  shell  can  be  made  by  comparing  the 

two-electron  electrostatic  integrals  FQ(ls,3p)  and  FQ(2p,3p) 

32 

obtainable  from  Mann's  tables  and  this  can  be  used  to 
modify  equation  1.43.  As  will  be  seen  later  AE^s  = 
(constant)  x  AE£p  to  a  very  good  degree  and  so  equation 
1.43  should  still  reflect  the  trend  of  the  Auger  shift  even 
if  the  potential  part  is  not  modified,  though  the  slope  may 
not  be  unity. 

The  transition  operator  gave  the  following  result  for 
the  two-step  process 


•  A  a  N , N— 1  n  N , N-2 . 

AE  (KLL  )  =  -A  (e  -  2e  ) 

rl 


(1.44) 


and  for  the  three-step  process 


24 


iEA(KLL'>  =  AE™  -  AeN-1'N"2 


(1.45) 


Application  of  the  "atom-charge  model"  to  equation  1.45  gave 


.  n  *  * 


■ 


' 


■ 
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kN-l,N-2 


The  charges  and  potentials  are  calculated  using  a  pseudo 


3 

atom  with  nuclear  charge  Z  +  j  in  place  of  the  original  Z 


atom.  The  parameters  being  obtained  from  interpolation  be¬ 
tween  those  associated  with  the  Z+l  atom  and  those  with  the 
Z+2  atom. 

The  three-step  model  is  the  easiest  to  handle  and  with 
the  findings  of  Adams  represents  the  more  obvious  application 
of  the  potential  model.  The  meaning  of  eN'N  ^  within  the 
TPM  formalism,  as  would  be  required  in  the  two-step  model, 
leads  to  possible  ambiguities  with  the  equivalent  core 
expression  for  the  singly-ionised  state  (sN  in  that  the 
same  parameters  for  the  central  atom  are  being  used  but 

33 

two  different  states  are  being  described.  It  was  noted 


that 


N,N-1  _ 


c  . 

l 


when  full  calculations  are  used  to  determine  the  orbital 
energies.  The  ambiguities  arise  when  the  potential  model 
is  applied.  These  can  be  removed  by  imposing  the  following 


relationships 


33 


«  , 
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A.N'N'2  =  lAto"'"-1  +  UN-1’N-2) 


where  oo  represents  the  "true"  potential.  By  assuming  a  sim¬ 
ilar  linearity  for  the  charges,  q,  potentials,  V,  and  K, 
namely 


N,N-1 

x 


N-l. 
x  ) 


N-l 


/  x 


1  ,  N,N-1 

2  x 


N-l , N-2 
+  x  ' 


) 


expressions  were  obtained  for  the  "correct"  equivalent  core 

potential^  (£) 


-  ztk*1'1*"1  -  kN)0 
N,N-2  =  ^N-l  +  (kN,N-l  _  kN)Q 


(1.47) 


,  „  N , N-l  N 

where  Q  =  q  -  q  . 

The  Auger  parameter,  a,  was  originally  defined  by 

tt  34 

Wagner  as 


a  =  (AEa  +  AEb)  =  2ARea 

with  AEg  being  a  core-binding  energy  shift  and  AR£A  the 

difference  in  extra  atomic  relaxation  of  the  initial  photo- 

33 

ionisation.  Application  of  equation  (1.47)  led  to 

A(Ea(KLL')  +  E  (L"))  =  2ARN'N"1  -  6AQ 

Arn,N  1  _  ^  ^~N  _  ^N'N  ,  the  molecular  relaxation  energy. 
This  expression  relates  the  shifts  to  the  ground  state  and 
transition  potential  of  the  initial  ionisation.  The  charge 


■  .  ■_  .V 


A. 


flow  term  (6AQ)  takes  account  of  the  differing  initial 
states  in  the  loss  of  the  first  electron  to  that  of  the 


second.  The  replacement  of  two-electron  term  by  terms 

associated  with  the  initial  photoionisation  removes  the 

ambiguities  however  it  is  necessary  to  infer  that  the 

relaxation  of  the  second  step  follows  from  the  first. 

33 

In  the  original  analysis  it  was  assumed  that 

Lr^K 

AE.,  was  close  to  zero.  This  assumption  can  lead  to 

substantial  errors.  The  ambiguity  can  be  avoided  by  using 
the  three-step  process  in  which  ^  is  used  in  con¬ 

junction  with  AE^_^y. 

A  more  accurate  expression  can  be  derived  for  the 
Auger  parameter. 


AEa  (KLL)  +  AEt  =  (AE  -AE  -AE  ,  )  +  AE  +  (AE  -AE  ) 
A  i-i  .Li  Uli 


L+K 

+  2AEt  -  AEt  - 

aet  , 

X-ray 

L  L 

L ' 

L+K 

X-ray 

-  A  (e  +e  )  + 

.  N-l ,N-2 

Ae 

.  N,N-1 

Ae  + 

,L-*K 

X-ray 

,  .  ,  N , N-l  Nv 

+  A  (e  -e  ) 

,  A  ,  N-l , N-2  N-l. 

+  A  (e  -e  ) 

,L-*K 

+  ARt  +  ARt  , 

(1.48) 

’X-ray 

L  L 

AEa(KLL)  +  AEk  =  (AE^ay- 

-ael.) 

+ 

AEk  +  (AEl-AEl) 

/  A 

AEx-ray 

-AEl,) 

+ 

<iEL+AEX-ray)  + 

(AEl-AEl) 

=  2AE^K 

+  art 

+ 

art  , 

(1.49) 
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where  the  orbital  approximations  for  the  binding  energy  have 
been  made  using  both  the  Liberman  and  transition  operator 
approaches.  The  final  expressions  indicate  that  the  shifts 
are  dependant  on  the  relaxation  accompanying  the  loss  of 
an  initial  L  shell  electron  (the  relaxation  accompanying 
the  loss  of  the  initial  K  electron  is  altered  upon  the 
level  being  filled  by  an  L  electron)  and  that  to  the  loss 

of  an  L  electron  leaving  an  ion  with  a  L-shell  hole.  Apply- 

..  .  .  —  —  .  n  ,  ,  .  N  N,N-1  N-l  ,  N-l ,N-2 

mg  the  kq  +  V  +  £  model  to  e  ,  e  ,  z  and  z 

will  give  estimates  of  RT ,  and  R  as  long  as  the  ap- 

propriate  values  for  k  are  known. 

The  Auger  parameter  formalism  provides  a  convenient 

way  of  relating  the  Auger  shift  to  the  binding  energy 

shift.  Rewriting  equation  1.48  gives 

T  ->-X 

AEAUG  =  (ARL  +  ARL->  -  (AEL  -  AEX-ray>  (1-50) 

This  is  nothing  more  than  a  rephrasing  of  equation  1.9.  It 

can  be  seen  that  the  Auger  chemical  shift  depends  on  a  ref- 

erence  level  (AET  -  AEV  )  which  is  modified  by  the 

L  A- ray 

relaxation  contribution  arising  from  the  loss  of  the  initial 
electron  and  that  of  the  second. 

D.  The  Purpose  of  This  Study 

As  with  core  electrons.  Auger  electrons  in  a  series  of 
molecules  also  show  a  chemical  shift.  Core  electron 
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binding:  shifts  have  been  studied  extensively,  however,  until 
recently,  there  has  been  very  little  systematic  study  of 
the  Auger  electron  shifts.  It  was  of  interest  to  study 
these  shifts  and  to  compare  them  to  the  photoelectron  shifts. 

To  have  readily  comparable  values  only  Auger  transitions 
involving  inner  core  electrons  can  be  considered.  The 
simplest  system  to  study  is  the  KLL  spectra  of  the  second 
row  elements.  The  complete  KLL  Auger  spectra  for  SF^,  S0o 

D  Z 

and  have  been  studied  quite  extensively . ^ ^ 

For  comparison  of  Auger  shifts  with  those  of  core 

electron  binding  energies,  only  the  KL^L^  line  need 

35 

be  considered.  Keski-Rahkonen  and  Krause  compared  the 

KL2^2  shift  with  the  Is  electron  shift  in  SF^,  H^S  and  SC^* 

3  6 

Asplund  et  al .  compared  the  ^2/2  with  the  Auger 

shift  in  SFg,  SC^ ,  ^S,  C0S  an<^  CS2*  In  ^is  latter  paper 

2  8 

a  formalism  based  upon  the  TPM  was  presented,  however  no 

comparison  of  theory  with  experimental  data  was  made. 

An  ab  initio  study  of  the  Auger  KL^L^  shift  compared 

.  .  31 

with  the  binding  energy  shifts  in  some  chlorine-containing 

37 

and  silicon-containing  species  was  conducted  by  Adams. 

18  2  5 

He  extended  the  Liberman  expression  '  for  binding 
energies  to  the  Auger  process  and  also  outlined  a  relaxation 


.8 


potential  expression. 

The  RPM  was  applied  by  Cavell  and  Sodhi^  to  a  series 
of  phosphorus  compounds  (PX^/  SPX^  /  OPX^  (X=F,C1);  PH^  and 
PF^)  for  both  the  core  electron  shifts  and  the  Auger  KL2L3 


. 
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shifts  using  a  CNDO/2  program.  The  approach  predicted 

the  Is  shifts  fairly  well  (and  to  a  lesser  extent  the  2p 

shift) .  Its  prediction  of  the  Auger  shifts  was  not  as 

satisfactory.  This  was  attributed  as  much  to  a  possible 

failure  of  the  CNDO  formalism  as  to  possible  limitations  in 

the  potential  model.  (Note,  after  publication  of  this  paper 

a  looping  error  in  the  part  of  the  program  calculating  the 

potential  at  the  nucleus  for  the  (Z+2)  equivalent  core  state 

was  discovered.  Thus  the  calculated  Auger  shifts,  as 

presented  in  the  paper,  are  in  error.  Corrected  shifts, 

as  appropriate  to  the  paper,  are  shown  in  appendix  3. 

The  conclusions  as  stated  above,  however,  are  still  valid.) 

The  major  feature  of  interest  was  the  lack  of  correlation 

between  the  Auger  shift  and  the  Pis  binding  energy  shift, 

not  only  were  the  magnitudes  different  but  in  some  cases 

there  was  a  reversal  in  trend.  This  was  also  observed  by 

3  6 

Asplund  et  al .  for  some  of  the  sulphur  compounds. 

24 

Kelfveetal.  have  performed  the  most  extensive 

previous  analysis  of  the  relationship  between  the  Auger 

shift  and  the  2p  electron  shift  for  a  series  of  silicon 

2  8 

compounds.  They  applied  the  TPM  using  a  modified  extended 

3  0 

Hiickel  program,  EWMO.  The  results  showed  satisfactory 
predictions  for  both  the  Auger  and  binding  energy  shifts. 

The  shifts  were  also  analysed  by  means  of  a  "group  shift" 
approach  and  in  the  limited  geometry  of  the  silicon  compounds 
gave  good  correlation. 


- 
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In  this  work  the  Is  and  2p  core  electron  shifts  are 
compared  to  the  KI^L^  (^D2^  Auger  electron  shifts  in  an 
extensive  series  of  phosphorus  and  sulphur  compounds.  Un¬ 
like  the  silicon  compounds  in  which  the  constant  geometry 
limits  general  conclusions,  compounds  of  sulphur  and  phos¬ 
phorus  show  various  coordination  environments.  It  was  of 
interest  therefore,  to  see  how  the  Auger  and  core  electron 
shifts  compared  and  hence  how  the  relaxation  effects  dif¬ 
fered  with  these  various  geometries.  By  studying  two  dif¬ 
ferent  series  of  compounds  it  is  possible  to  see  how  the 
effects  of  one  substituent  transfer  from  one  centre  to 
another.  We  have  therefore  tested  various  phosphorus  and 
sulphur  compounds  and  analysed  the  results  in  terms  of  the 
various  models  available. 


Chapter  2.  Experimental 

The  Is  and  2p  core  levels  and  the  KL2L3  Auger  line 
were  recorded  using  a  McPherson  ESCA  36  photoelectron 
spectrometer. 

The  2p  line  was  produced  by  irradiation  with  the  Al 

K  (1486.65  eV)  X-ray  or  the  Mg  K  (1253.64  eV)  X-ray  and 
ut  cx 

referenced  to  the  Kr  3p^  /0  and  Kr  3(1^.  /n  lines  (214.55  (15) 


3/2 


‘5/2 
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eV  and  93.80  (10)  eV  binding  energies  respectively). 

The  Is  line  was  produced  by  irradiation  with  the  Ag 

35 

(2984.34  (2)  eV)  X-ray  except  for  oxygen  containing 
sulphur  compounds,  in  which  case  the  Ag  LD ,  (3150.97  (3J  eV) 

p  i 

X-ray  was  used.  The  KI^L^  line  was  collected  in  the  same 

scan  as  the  Is  line  and  both  lines  were  referenced  to  the 

4 1 

Ne  KL2  3K2  3  Auger  line  (804.56  (2)  eV  kinetic  energy)  and 
the  Ne  Is  photoelectron  line  (870.37  (9)  eV  binding  energy) ^ 
produced  by  the  Ag  L  ^  X-ray. 

In  all  cases  the  X-ray  sources  were  operated  at  a 
power  of  400  watts  (10  kV  x  40  mA) . 

The  reference  gas  supply  was  controlled  with  a 
Granville -Phi Hips  series  203  variable  leak  valve.  Liquid 
samples  were  contained  in  an  evacuated  glass  vial  equipped 
with  a  Teflon  valve.  The  vapour  was  introduced  via  a 
Nupro  SS-4BMW  all  metal  bellows  metering  valve.  Gaseous 
samples  contained  in  cylinders  were  introduced  via  either 
the  Granville -Phi llips  or  the  Nupro  valve.  In  order  to 


42 


, 


■ 


43 


control  the  gas  pressure  accurately  the  cylinder  pressure 
was  reduced  by  a  Matheson  3500  series  regulator. 

All  lines  were  made  of  1/4"  stainless  steel  tubing 
and  connections  were  made  with  1/4"  Swagelok  fittings.  The 
sample  (and  reference)  supplies  were  mixed  in  the  feed 
line  prior  to  introduction  into  the  sample  cell.  X-rays 
entered  the  sample  cell  through  an  aluminum  window 
(0.0001"  thickness). 

The  pressure  of  the  mixture  was  monitored  using  a 
MKS  Baratron  Pressure  meter.  The  sample  and  reference 
gases  were  of  approximately  equal  pressures  and  typical 
(total)  pressures  ranged  from  150-200  u.  The  mixtures  were 
successively  scanned  until  adequate  statistics  (1500  -  c/s 
on  highest  peak)  were  achieved.  Repeated  rapid  cycles  were 
run  to  minimise  any  possible  error  between  reference  and 
sample  peaks  due  to  slight  pressure  drift. 

The  kinetic  energy  of  the  electron  is  calculated  from 
the  equation 


K.E.  =  k  AV 

where  AV  is  the  potential  difference  between  the  analyser 
plates  and  k  is  the  experimentally  determined  machine  constant. 
In  order  to  compensate  for  any  variations  with  time  in  the 
machine  constant,  k,  which  are  noticable  over  the  large 
kinetic  energy  differences  involved  in  this  work,  all  peak 
positions  were  normalised  to  give  reference-line  separations 


. 

- 


in  agreement  with  those  given  in  the  literature,  ie . 


[reference  separation]  [x  -  reference] 

cAU  L  ___  cXp  t 


[reference  separation] ^ 


[x  .  ,  -  reference] . 

actual  lit 


hence  the  true  kinetic  energy  of  x,  xactuai'  was  found. 
Subtracting  the  core  electron  kinetic  energies  from  the 
energy  of  the  exciting  source  gave  the  required  binding 
energy. 

At  least  three  separate  runs  were  performed  for 
each  data  set  and  an  average  value  was  taken.  The  largest 
deviation  from  the  average  was  0.09  eV  with  the  majority 
of  cases  deviating  by  0.05  eV  or  less.  It  is  concluded 
that  the  errors  in  the  shifts  are  no  more  than  ±0.1  eV. 

Positions  and  areas  of  each  peak  were  obtained  by 

42 

using  a  non-linear  least  squares  program  which  could  fit 
either  a  Lorentzian  or  Gaussian  curve.  No  attempt  was 
made  to  deconvolute  the  spin  orbital  components  of  2p  peak 
as  the  resolution  was  insufficient.  Some  difficulty  was 
experienced  in  deconvoluting  the  Kr  3d  peak  however  con¬ 
straining  the  Kr  an<3  Kr  ^3/2  to  e<3ual  gave  a 

good  fit  with  the  required  intensity  ratios  and  separations 
A  sufficient  portion  of  the  low  energy  side  of  the  KI^L^ 
Auger  line  was  collected  in  order  to  ensure  a  good  fit  on 
the  KL2L3  Auger  peak,  however  this  portion  of  the  spectrum 
was  not  accurately  fitted. 
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All  the  compounds  were  obtained  commercially  with  the 

43  44 

exception  of  SPF^  and  PCCF^)^  which  were  made  by 

standard  literature  procedures  in  this  laboratory,  and 

P(C2H^)^  which  was  supplied  by  Dr.  John  Malito  of  this 

department.  Where  necessary,  the  reagents  were  fractionated 

under  vacuum  to  remove  impurities,  and  the  purity  of  the 

compound  was  verified  by  I.R.  or  N.M.R.  spectroscopy  or 

45 

both.  Chemical  means  had  to  be  employed  to  purify  SF^  , 
otherwise  the  compounds  were  run  without  further  purification. 
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Chapter  3.  Calculations 


A.  Semi-empirical  Programs 

Two  types  of  semi-empirical  calculations  were  em- 

•50 

ployed.  The  first  utilised  the  CNDO/2  formalisation 

using  the  original  parametrisation  of  Pople  et  al.  with 

the  exception  of  Cl.  The  second  was  a  modified  iterative 

extended  Hiickel  program  (EWMO)  which  was  adapted  for 

transition  potential  calculations.  Calculations,  in  both 

cases,  were  performed  with  d-orbitals  not  included  in 

the  basis  set.  In  the  CNDO/2  case  attempts  were  made  to 

include  d-orbitals,  however  difficulties  occurred  with 

convergence  for  many  of  the  molecules. 

Both  programs  had  to  be  parametrised  for  Ar  so  that 

calculations  could  be  performed  to  estimate  the  Auger 

electron  chemical  shift  for  the  sulphur  compounds. 

46 

In  the  CNDO/2  program  the  new  Cl  parameters  '  were 

1/2 (I  +  A)  =  -21.76  ev 
s 

1/2 (I  +  A)  =  -  9.56  eV 

The  Ar  parameters  used  for  the  CNDO  calculations  were 

,  47 

those  estimated  by  St/gard  and  Manne . 

The  EWMO  (Hiickel)  program  has  been  used  successfully 
by  Kelfve  et  al. 24  on  a  series  of  silicon  compounds.  The 
main  features  of  the  program  have  been  adequately  described^ ' 4^ 
and  will  not  be  dealt  with  here. 
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The  Ar  parameters  for  the  EWMO  (Hiickel)  program  were 

obtained  in  the  following  manner.  The  expression  for  the 

electron-repulsion  integrals,  g  ,  within  the  CNDO  form- 

s  s 

alisation  reduces  to 


ss 


=  [SASA 


S  K  ] 
A  A 


where  represents  the  orbitals  and  are  Slater-type  s 
orbitals.  This  reduces  to 


g  =  const,  l 
^ss 

where  £  is  the  Slater  exponent.  The  ratio,  g  (Ar) /g  (Cl) 

s  s  s  s 

then  becomes  the  ratio  of  the  Slater  exponents: 

g  (Ar)  c  2.25 

-  =  =  _  =  1  11 

g  (Cl)  «C1  2.03  ‘ 

By  using  the  values  of  the  average  one-centre  electron 

repulsion  integrals  anc^  the  average  two-centre 

* 

electron  repulsion  integrals  (g  )  appropriate  to  the 

4  9 

EWMO  (Hiickel)  calculations  for  Cl  ,  the  following  estimates 
for  Ar  were  obtained: 

g,  _  =  11.42  eV  g*  =  11.51  eV 
^AA  AA 


These  values  and  the  experimental  ionisation  potentials 
for  Ar  3s  and  Ar  3p  (29.24  eV  and  15.82  eV  respectively) 


give  the  core-matrix  elements,  Usg  and  upp: 


"  •'  '  . 


- 
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U  =  -109.78  eV  U  -96.36  eV 

ss  pp 

by  use  of  the  relationship 


U.  .  + 

li 


<aA-1J 


In  all  cases  calculations  were  performed  on  each 
molecule  for  the  ground  state  (Z)  and  the  equivalent  core 
states  (Z+l)  and  (Z+2) .  The  equivalent  core  states  were 
obtained  by  replacing  the  atom  of  interest  with  the  (Z+n) 
atom.  A  molecular  charge  of  n  throughout  the  series 
maintained  a  constant  number  of  electrons.  In  this  way 
the  properties  of  the  singly-ionised  and  doubly-ionised  ions 
were  mimicked. 

Furthermore,  in  the  Hiickel  case,  calculations  were 
performed  on  each  molecule  for  the  transition  states 
(Z+l/2)  and  (Z+3/2).  In  this  case  the  atom  of  interest 
was  replaced  by  a  pseudo-atom.  For  the  (Z+l/2)  case  the 
pseudo-atom  parameters  were  obtained  by  interpolation 
between  those  of  the  Z  atom  and  those  of  the  (Z+l)  atom. 

In  the  (Z+3/2)  case  the  parameters  were  interpolated  be¬ 
tween  the  (Z+l)  atom  and  the  (Z+2)  atom. 


B.  Models  for  Chemical  Shift 

Semi-empirical  calculations  on  ESCA  shifts  have  proved 
useful  in  the  context  of  the  potential  model.  The  original 
form  of  the  potential  model  is 


-• 
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AE_  =  k.qA  +  V.  +  £ 
B  i^A  A 


(3.1) 


This  will  be  termed  the  "atom  charge  model". 

The  terms  in  equation  (3.1)  are: 

AEB  is  the  binding  energy  shift. 

is  the  charge  on  atom  A,  the  atom  of  interest. 

V  =  E  qn/R^  ,  the  potential  on  the  ionised  atom,  A,  due 
A  B^A  B  AB 


to  the  charges  at  all  other  atoms  in  the 
molecule . 

k.  =  <ij | i/r .  .  | ij>  ,  the  coulomb  repulsion  integral  between 
l  13 

the  core  electron  (i)  and  the  valence 
electron  ( j ) . 

£  is  a  constant  depending  on  the  reference  level. 

In  practice  and  £  are  obtained  by  a  least-squares 
fit  on  experimental  data.  This  allows  a  value  of  k  for 
each  core  level  to  be  estimated.  However,  in  practice, 
it  is  found  that  the  value  of  k  varies  according  to  the 
compounds  selected  for  the  least-squares  fit. 

By  invoking  the  point-charge  approximation,  k.  becomes: 


=  <j|i/ril j> 


£ 

n 


(3.2) 


that  is  the  nuclear  attraction  integral  of  the  valence 
orbital  at  A  for  nucleus  A,  which  in  turn  is  equal  to  the 
Slater  orbital  exponent  (c)  for  the  valence  shell  divided 
by  the  principal  quantum  number  of  the  valence  shell. 


. 
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Equation  (3.1)  can  be  modified  to  form  the  potential 
model : 


AEb  =  -AV(Z)  =  -A(KPaa  -  ^  qB/RAB)  (3.3) 

where  is  the  valence  electron  population  on  atom  A  and 

K  is  calculated  for  each  case  using  the  point  charge  ap- 

23a 

proximation.  This  is  the  approach  used  by  Davis  et  al. 

Use  of  the  ground  state  (i.e.  parameters  associated 
with  the  Z  atom)  in  equation  (3.1)  or  (3.3)  allows  an 
estimation  of  the  chemical  shift  to  be  obtained.  This 
approach  ignores  the  relaxation  occurring  upon  photo¬ 
ionisation  and  is  known  as  the  Ground  Potential  Model  (GPM) . 

The  Relaxation  Potential  Model  (RPM)  attempts  to 
correct  for  relaxation  by  use  of  the  equivalent  core 
approach.  Equation  (3.1)  is  modified  by  using  charges  and 
the  off-centre  atom  potentials  given  by 


qA  =  l/2(qA(Z)  +  qA(Z+l)) 
VA  =  1/2 (VA(Z)  +  VA(Z+1) ) 


(3.4) 


The  appropriate  k  is  obtained  theoretically  or  by  fitting. 

23a 

The  RPM  as  used  by  Davis  etal.  is  obtained  by 
taking  the  negative  means  of  V(Z)  and  V(Z+1)  found  by  use 
of  equation  (3.3).  That  is 


AEd  =  A(V(Z)  +  V(Z+1)  ) 

D  2. 


(3.5) 
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The  use  of  the  pseudo-atom  (the  Transition  Potential 
Model  -  TPM)  allows  a  direct  estimate  of  the  relaxed 
binding  energies  from  equation  (3.1).  The  appropriate 
charges  and  off-atom  potentials,  in  this  case,  are  given 
in  one  calculation. 

An  estimation  for  the  relaxation  energy  can  be  ob¬ 
tained  by  subtracting  the  GPM  value  for  the  chemical  shift 
from  the  RPM  or  TPM  value. 

The  extension  into  Auger  spectroscopy  is  obtained  by 

using 

AEAUG  =  A(EB  (K)  -  Eb  (L)  -  Eb  (L*  )  )  =  AExray  -  AEb  (L*  )  (3.6) 


X 

where  AEAUG  is  the  KLL  Auger  shift.  The  term  AEg (L  )  is 
the  binding  energy  shift  for  an  L-shell  electron  in  a 
species  already  containing  an  L-shell  hole. 

The  appropriate  models,  as  described  earlier,  can  be 
used  by  utilising  a  (Z+l)  initial  state  and  a  ( Z+2)  final 
state  or  a  (Z+3/2)  transition  state. 

In  the  CNDO  calculations  ground  and  relaxation  po¬ 
tential  models  were  utilised  for  both  the  atom-charge  and 
potential  models.  In  the  former  case  both  fitted  and 
theoretical  k's  were  used.  In  the  Hiickel  case  the  ground, 
relaxation  and  transition  potential  models  were  utilised 


for  the  atom-charge  model. 

All  the  geometries,  where  obtainable,  were  experimental. 


■u  '  “i  '  1  3 


Other  geometries  were  assumed  using  bond  lengths  and 
angles  from  similar  compounds  of  known  geometry.  Full 
details  of  the  geometries  used  in  the  calculations  can 
be  found  in  Appendix  1. 


- 
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Chapter  4.  Results  and  Discussion 


A.  Experimental  Results 

All  shifts  were  referenced  to  PH^  in  the  phosphorus 
series  and  to  H^S  in  the  sulphur  series.  Phosphorus  Is  and 
2p  photoelectrons  and  KL^L^  Auger  shifts  are  given  in 
Table  4.1  and  the  corresponding  data  for  the  sulphur  series 
in  Table  4.2.  The  correlation  between  the  Is  and  2p 
photoelectron  shifts  for  phosphorus  is  shown  in  Figure  4.1 
and  the  correlation  between  the  Is  photoelectron  and  KI^L^ 
Auger  shifts  for  phosphorus  is  shown  in  Figure  4.2.  The 
corresponding  correlations  for  sulphur  are  shown  in  Figures 
4 . 3  and  4.4. 

It  was  found  that  the  oxygen  KW  Auger  region  coincided 
with  the  Sis  line  under  Ag  L  ^  X-ray  excitation.  Ac¬ 
cordingly  Sis  peaks  were  collected  with  Ag  X-ray 

excitation  for  the  oxygen-containing  sulphur  compounds.  The 
Sis  binding  energies  for  H^S  differed  when  Ag  X-ray 

or  Ag  LD ..  X-ray  excitation  were  compared;  the  former  gave 
pi 

a  Sis  value  of  2478.87  eV  whereas  the  latter  gave  a  value 
of  2478.60  eV.  Since  the  Ag  Lal  X-ray  gave  a  more  intense 
photoelectron  line  with  a  smaller  FWHM  than  Ag  X-ray 

the  former  value  was  considered  to  be  more  reliable.  All 
Sis  peaks  obtained  by  means  of  Ag  X-ray  excitation  were 

referenced  to  the  Sis  line  of  produced  by  the  Ag 

X-ray  in  order  to  obviate  binding-energy  shift  errors. 
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Table  4.1. 

Experimental 

chemical 

shifts3 

relative 

to  PH3b 

Compound  # 

AlsC 

A2p^ 

-akl2l3 

PH3 

1 

0.00 

0.00 

0.00 

spf3 

2 

6.20 

5.15 

2.98 

p(ch3)3 

3 

-0.98 

-1.10 

-1.0S® 

-1.11 

-3.30 

spci3 

4 

4.21 

3.57 

52 

3 . 58DZ 

-0.47 

opf3 

5 

6.96 

5.91 

52 

6.03  z 

5.14 

0PC13 

6 

4.70 

4.01 

52 

4 . 01DZ 

0.72 

PF5 

7 

8.55 

7.33 

52 

8.83 

5.93 

PC13 

8 

3.28 

2.79 

2.7352 

-0.39 

PF3 

9 

5.48 

4.72 

4.62 

4.76 

SP(CH30)  Cl 

10 

3.60 

2.78 

-0.95 

(ch3o)3po 

11 

3.34 

2.54 

-0.11 

(CHjO) 3PS 

12 

3.26 

2.39 

-0.95 

(ch30)3p 

13 

1.82 

1.19 

-1.03 

P (CH3)C12 

14 

1.89 

1.53 

-1.33 

P (C2H5  >  3 

15 

-1.48 

-1.53 

-4.48 

OP  (CH  C1)C1 

2  16 

3.77 

3.18 

-0.10 

sp(ch3)ci2 

17 

3.19 

2.59 

-1.44 

p(cf3)3 

18 

1.74 

1.70 

-1.16 

(a)  All  shifts  taken  as 

Ashift  = 

compound 

-  reference. 

Shifts 

are  reliable 

to  ±0.1  < 

sV. 

(Cont ' d . ) 


, 
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Table  4.1.  Cont'd. 


(b)  Reference  lines  for  PH^  (in  eV) 


This  work 


Ref.  51 


Ref.  52 


Is 

2p 

KL2L3 


2150.88  (20) 
137.35  (20) 
1841.46  (20) 


137.3 

1841.4 


136.87 


(c)  The  Is  lines  were  produced  by  AgL  (2984.34  eV) 

0 it 

Auger  line  was  collected  in  the  same  scan. 


Calibration  lines  used  were: 


Ne  KL23L23  804.56  (2)  eV 

AgL^  converted  by  Ne  Is  870.37  (9)  eV 

(d)  The  2p  lines  were  produced  by  AIK  (1486.65  eV) 

Uw 

MgKa  (1253.64  eV) . 

Calibration  lines  used  were: 


Kr  3d5/2 

93.80 

(10)  eV 

Kr  3p3/2 

214.55 

(15)  eV 

(e)  Values  adjusted  from  Ref.  51. 


.  The 


or 
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Table  4.2.  Experimental  chemical  shifts  relative  to 
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Figure  4.4  Experimental  (-)KL2L3  Auger  energy  shifts  versus  Is  binding  energy 
shifts  -  sulphur  series.  Compound  numbers  are  presented  in  Table  4.2. 
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Agreement  between  the  shifts  reported  here  and  those 
of  other  workers  was  found  to  be  good  with  the  exception 

52 

of  some  of  the  2p  shifts  shifts  reported  by  Jolly  et  al. 

The  largest  discrepancy  was  1.5  eV  for  the  P2p  shift  of 

PF,_ .  Our  PFj.  2p  spectrum  was  redetermined  on  different 

occasions  and  we  are  confident  of  the  present  value.  Further 

57 

support  for  the  view  that  the  value  quoted  by  Jolly  et  al. 

is  in  error  is  that  their  value  gives  a  larger  PH^  to 

PF,-  shift  for  the  2p  shift  than  that  for  our  Is  shift  which 

is  inconsistent  with  the  general  ls/2p  shift  relationship. 

The  absolute  values  of  the  photoelectron  and  Auger 

lines  for  the  PH^  reference  are  given  in  Table  4.1.  The 

50 

value  of  Ashe  et  al.  for  the  P2p  line  (137.3  eV)  agrees 
well  with  the  value  of  137.35  eV  found  in  this  work,  whereas 
the  value  quoted  by  Jolly  et  al.  (136.87  eV)  is  sig¬ 
nificantly  lower.  The  KL^L^  auger  value  of  Ashe  et  al.50 
(1841.4  eV)  also  agrees  well  with  our  value  of  1841.46  eV. 

The  absolute  values  for  the  H2S  reference  are  given 
in  Table  4.2.  Agreement  with  the  S  KL2L3  Auger  value  of 
Asplund  et  al.  is  excellent  as  is  the  agreement  for  the 
S2p  photoelectron  value,  once  the  spin-orbit  coupling  is 
taken  into  account.  A  systematic  calibration  error 
(v.4  eV)  seems  to  be  present  in  the  Sis  photoelectron 
energy  between  this  work  and  that  of  Keski-Rahkonen  and 
Krause,  however,  the  reported  shifts  agree  well  within 
the  limits  of  experimental  error. 
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As  can  be  seen  in  Figures  4.1  and  4.3,  the  cor¬ 
relation  between  the  Is  and  2p  photoelectron  shifts  is  good 
with  the  Is  shift  being  larger,  however,  when  comparing  the 
photoelectron  with  the  Auger  shifts  (Figs.  4.2  and  4.4), 
there  is  poor  correlation.  The  photoelectron  and  Auger 
shifts  of  an  atom  in  a  series  of  compounds  of  the  same 
coordination  state  would  be  the  most  likely  series  to  show 
any  correlation  of  these  values.  Comparisons  of  compounds 
in  which  the  coordination  number  of  the  central  atom  is 
altered  is  the  least  favourable  case  for  correlation. 

While  the  latter  point  is  supported  by  the  data  (compare 
the  relationship  between  S02F2  and  SF^Cl  or  SPF^  and  PF^) 
similar  differences  also  occur  in  comparisons  of  the  former 
type,  the  Auger  shift  between  SO(CH2)2  and  S0C12  is  smaller 
than  the  corresponding  Sis  shift  whereas  the  converse  is 
true  in  the  case  of  SOCl2  and  SOF2 .  Another  example  is  in 
the  shifts  of  P(CF2)^  and  PCI,  from  PH^ .  The  differences 
in  trends  must  arise  from  the  relaxation  contribution 
associated  with  the  initial  photoionisation  of  the  molecule 
and  that  associated  with  the  removal  of  a  second  electron 
from  the  ion.  It  is  this  point  that  the  various  theoretical 
models  must  resolve. 

B.  Application  of  the  Experimental  Results  to  the  Group 

Shift  Concept 


Group  shift  concepts  might  provide  a  useful  means  of 
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predicting  chemical  shifts  in  Auger  and  photoelectron 
spectroscopy.  It  involves  the  partioning  of  the  total 
shift  into  contributions  made  by  each  group.  Inherent 
in  the  group  shift  idea  is  the  assumption  that  shifts 
incurred  by  each  group  are  independent  of  each  other. 

This  assumption  ignores  the  fact  that  more  electronegative 
substituents  can  induce  a  higher  charge  on  the  central 
atom  which  in  turn  can  affect  the  other  substituents. 

A  test  of  the  group  shift  concept  on  a  series  of 
silicon  compounds  indicated  that  there  was  a  very  good 

24 

correlation  between  group  shift  and  total  chemical  shift. 

In  the  silicon  series  however,  all  the  molecules  encountered 
were  tetrahedral  and  only  a  limited  range  of  groups,  -Cl, 
OCI^CH^,  -H,  -CH^ ,  -CH=CH2  and  -CI^CH^  were  evaluated.  In 
the  phosphorus  and  sulphur  compounds  presented  here,  ad¬ 
ditional  factors  such  as  the  change  of  the  coordination 
number,  the  effect  of  lone  pairs,  the  inclusion  of  highly 
electronegative  groups  such  as  -F  and  -CF^  and  groups 
bonded  to  the  central  atom  with  a  double  bond  were  in¬ 
troduced  which  should  provide  a  more  thorough  test  upon 
the  general  applicability  of  the  group  shift  model. 

Experimental  group  shifts  (group  x  -  group  y)  are 
given  in  Table  4.3.  It  is  apparent  that  in  spite  of 
occasional  good  agreement  there  exists  serious  discrepancies. 

For  the  Pis  series  the  agreement  is  good  with  the  exception  of 
any  compounds  involving  a  terminal  sulphur  substituent. 
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For  the  P  KL^L^  series  there  are  many  more  exceptions  and 
the  validity  of  the  concept  appears  to  have  vanished.  Thi 
view  is  further  supported  by  the  values  for  the  series  of 
compounds  with  sulphur  as  the  central  atom.  Comparison  of 
group  shifts  obtained  from  compounds  of  differing  co¬ 
ordination  indicate  the  lack  of  general  applicability 
readily  revealed  by  comparing  the  shift  obtained  from  the 
difference  ((CH^)  2^0- (CH^)  2^)  to  that  obtained  from  the 
difference  (  (CH^O)  2^2“  (CH^O)  2^0)  .  Furthermore  the 
(CH^-H)  difference  illustrates  the  lack  of  consistency  in 
the  group  shift  method  even  when  compounds  of  the  same 
coordination  number  are  compared. 

To  summarise,  the  group  shift  concept  appears  to  be 
valid  only  within  very  limited  series  of  compounds  at  best 
and  it  must  be  used  with  caution,  hence  the  concept  is  of 
little  use. 

C.  Results  of  Calculations  -  Charges  and  Potentials 

All  of  the  calculations  performed  gave  the  net  charge 
on  the  central  atom,  q^,  and  the  potential  on  the  central 
atom  due  to  the  charges  at  all  other  atoms  in  the  molecule 
VA.  Charges  and  potentials  were  obtained  for  the  neutral 
molecule  (q  N,  V  N) ,  its  singly-ionised  core-hole  state 

v  N_1)  and  its  doubly-ionised  core-hole  state 
A  A 

(q  N“2,  v  N“2).  To  obtain  values  for  the  ionised  species 
A  A 
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within -the  context  of  the  semi-empirical  calculations  used 

here  the  equivalent-core  approximation  was  used.  The 

central  (Z)  atom  in  the  molecule  was  replaced  by  the  iso- 

electronic  (Z+l)+  and  ( Z+2 )++ cations  for  the  singly-ionised 

and  doubly-ionised  species  respectively.  Transition  state 

charges  and  potentials  for  a  N  -*■  (N-l)  and  a  (N-l)  (N-2) 

transition  (qA  '  Va  ;  qA  '  VA  )  were 

also  obtained,  but  only  with  the  EWMO  (Hiickel)  program. 

+1/2 

The  central  (Z)  atom  was  replaced  by  a  (Z+l/2)  pseudo-atom 


for  the  N  (N-l)  transition  and  by  a  (  (Z+l)  +1/2) 


+  3/2 


pseudo¬ 


atom  for  the  (N-l)  -*  (N-2)  transition.  Parameters  for  the 
pseudo-atom  were  obtained  by  interpolation  between  the  Z, 
Z+l  and  Z+2  atomic  parameters. 

Results  of  the  CNDO/2  calculation  are  given  in  Tables 
4.4  and  4.5  for  the  phosphorus  and  sulphur  series  res¬ 
pectively.  The  corresponding  results  of  the  EWMO  (Hiickel) 
calculations  are  given  in  Tables  4.6  and  4.7.  Convergence 
problems  occurred  in  the  EWMO  program  in  the  case  of 
P(CF3)3  and  for  the  (N-2)  state  of  P(C2H5)3  thus  these 
values  are  not  available.  In  all  other  cases  convergence 
was  achieved. 

The  net  charge  on  the  atom  being  ionised  (qA)  is 
related  to  its  valence  population  (P^)  by 
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Table  4.4.  Central  Atom  Charge  (q)  and  Off-atom  Potential  (V)  From  CNDO/2  Cal- 
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Table  4.5.  Central  atom  charge  (q)  and  off-atom  potential  (V)  from  CNDO/2  Calculation 
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Table  4.6.  Central  atom  charge  (q)  and  off-atom  potential  (V)  from  EWMO  (original  parameterisation) - 
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Table  4.7.  Central  atom  charge  (q)  and  of f-atom  potential  (V)  fromEWMO  (original  parameterisation)  -  S  series 
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* 

where  Z^  is  the  effective  reduced  charge  of  atom  A  (the 
nuclear  charge  minus  the  number  of  core-electrons) .  Plots 
of  the  phosphorus  valence  population  versus  the  state  of 
core-ionisation  derived  from  CNDO/2  calculations  are  shown 
in  Figure  4.5  and  from  EWMO  (Hiickel)  calculations  in 
Figure  4.6.  Similar  plots  for  the  sulphur  series  are 
shown  in  Figures  4.7  and  4.8.  The  change  in  valence  pop¬ 
ulation  upon  core-ionisation  gives  a  measure  of  the  electron 
flow  towards  the  central  atom  due  to  the  formation  of  a 
hole  state. 

(i)  CNDO/2  Results 

The  CNDO/2  calculation  gives  an  almost  linear  increase 
in  the  valence  population  of  the  core  atom  upon  increasing 
core-ionisation.  A  larger  increase  in  the  valence  population 
on  the  core-atom  is  perceivable  upon  second  ionisation  for 
some  of  the  molecules,  however  for  all  practical  purposes 
the  trend  is  linear.  A  similar  linear  trend  was  observed 

O  A 

by  Kelfveetal.  for  CNDO/2  calculations  performed  on  the 
series  SiCl  (CH~)„  ,  however  the  calculation  predicted 

that  SiCl^  would  have  the  highest  silicon  valence  pop¬ 
ulation  in  the  series,  contradicting  the  results  expected 

24 

on  the  basis  of  simple  electronegativity  arguments. 

This  behaviour  was  not  observed  in  the  present  series  of 
calculations  for  phosphorus  or  sulphur  compounds.  The 


. 


... 

;  ■ 


j  1 


. 


. 


77 


Figure  4 . 5  Comparison  of  the  P  valence  population  (from 
CNDO/2  calculation)  with  the  state  of  core  ionisation  for 
some  phosphorus  compounds.  Compound  numbers  are  presented 


in  Table  4.1. 
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State  of  Core  Ionisation 

Figure  4 . 6  Comparison  of  the  P  valence  population  (from 
EWMO  calculation)  with  the  state  of  core  ionisation  for  some 
Phosphorus  compounds.  Compound  numbers  are  presented  in 


Table  4.1. 


.  ' 


■ 


79 


G 

O 

•H 

-P 

td 

i—i 

3 

a 

O 

Ph 

(D 

U 

G 

03 

i— I 

tO 

> 


Figure  4,7  Comparison  of  the  S  valence  population  (from 
CNDO/2  calculation)  with  the  state  of  core  ionisation  for 
some  sulphur  compounds.  Compound  numbers  are  presented  in 


Table  4.2. 
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Figure  4 . 8  Comparison  of  the  S  valence  populatioi 


(from 


EWMO  calculation)  with  the  state  of  core  ionisation  for  some 
sulphur  compounds.  Compound  numbers  are  presented  in  Table  4.2. 


' 


valence  population  of  the  core  atom  in  both  cases  follows 
the  expected  electronegativity  trends. 

The  CNDO/2  calculations  performed  by  Kelfve  et  al. 
included  d-orbitals  whereas  the  calculations  reported 
here  did  not.  For  comparison  CNDO/2  calculations  including 
d-orbitals  were  performed  on  PCl^  and  PCCH^)^.  The  results 
are  presented  in  Table  4.8  and  Figure  4.9.  The  expected 
trend  in  the  phosphorus  valence  population  is  preserved 
when  comparing  the  population  in  PCCH^)^  to  PCl^,  however 
the  calculation  shows  more  of  an  increase  in  the  valence 
population  on  the  core  atom  upon  second  ionisation  compared 
to  the  initial  ionisation  than  in  the  case  where  d-orbitals 
were  omitted  from  the  basis  set.  These  changes  arising 
from  a  change  in  basis  set  illustrate  the  dangers  of  ob¬ 
taining  charges  and  potentials  from  semi-empirical  methods 
involving  a  minimal  basis  set. 

(ii)  EWMO  (Hlickel)  Results 

The  EWMO  (Huckel)  calculation  shows  different  trends 
for  the  change  in  valence  population  of  the  core  atom 
upon  increasing  core  ionisation  in  the  phosphorus  and 
sulphur  series.  The  sulphur  series  shows  an  almost  linear 
relationship  of  population  with  possibly  an  indication  of 
saturation.  The  phosphorus  series  shows  a  markedly  larger 
increase  in  the  valence  population  in  going  from  the  first 
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Table  4.8.  Phosphorus  valence  population  on  P(CH^)^  and 

PCl^  from  CNDO/2  calculations  with  d-orbitals 
included . 

State  of  Core  Ionisation 


N 

N— 1 

N-2 

PC13 

4.7338 

5.4891 

6.4732 

P<CH3>3 

4.9113 

5.6569 

6.5815 
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Figure  4.9  Comparison  of  the  P  valence  population  (from 
CNDO/2  calculation)  with  the  state  of  core  ionisation  for 
PCCH^)^  an^  PCI3  with  and  without  d  orbitals  included  in 
the  basis  set. 

G  PCl^  with  d  orbitals  in  basis  set 
A  PCl^  without  d  orbitals  in  basis  set 
S  P (01^)3  with  d  orbitals  in  basis  set 
O  P (CH^ ) ^  without  d  orbitals  in  basis  set 


, 
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ionisation  to  the  second.  Both  of  these  trends  contrast 

O  A 

with  that  shown  by  Kelfveetal.  for  silicon.  Upon  in¬ 
creasing  core-ionisation,  the  silicon  series  showed  a 
marked  saturation  in  the  increase  of  the  valence  pop¬ 
ulation  on  the  central  atom. 

This  difference  in  trends  between  the  three  series 
of  compounds  for  the  same  type  of  calculation  is  disturbing. 
A  possible  explanation  for  this  occurrence  might  be  in  the 
parameterisation  of  the  phosphorus  atom  in  the  EWMO  (Hiickel) 
program.  As  a  test  to  see  whether  the  behaviour  of  the 
phosphorus  compounds  is  a  quirk  of  the  parameterisation  or 
whether  it  is  part  of  a  general  trend,  a  few  of  the 
phosphorus  compounds  were  run  with  the  silicon  anion  (Si  ) 
as  an  equivalent  core.  The  results  are  shown  in  Figure 

4.10.  It  can  be  seen  that  the  values  involving  the  use 

of  the  phosphorus  parameters  are  too  high  and  this  implies 
that  the  trend  in  the  phosphorus  series  is  caused  by  the 
particular  parameterisation  chosen  for  phosphorus  in  the 
EWMO  (Hiickel)  program. 

In  view  of  this  result  it  was  decided  to  change  the 
parameters  for  phosphorus  in  the  EWMO  (Hiickel)  program. 

Full  details  of  the  methods  used  are  given  in  Appendix  2. 

A  plot  of  SiCl^  valence  population  versus  the  state  of 
core-ionisation  obtained  with  both  the  original  para¬ 
meterisation4  8  and  for  the  new  parameters  is  given  in  Figure 

4.11.  For  comparison,  similar  data  for  PCl^  and  SC^  are 
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State  of  Core  Ionization 


Figure  4.10  Comparison  of  the  P  valence  population  (from 
EWMO  calculation)  with  the  state  of  core  ionisation  for  some 
phosphorus  compounds  going  from  the  Si  equivalent  core 
state  to  the  Cl++  equivalent  core  state. 

#  opf3  O  pf3  □  opci3  ©  pci3 
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Figure  4.11  Comparison  of  the  central  atom  valence  pop¬ 
ulation  (from  EWMO  calculation)  with  the  state  of  core 

ionisation  for  SiCl.,  PC1_  and  SC10. 

4  3  2 

O  original  P  parameterisation  <^>  new  P  parameter- 

isation 
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also  given. 

It  can  be  seen  that  all  three  series  of  compounds 
now  proceed  in  a  similar  manner  upon  increasing  core¬ 
ionisation.  Also,  the  transition  potential  model  values 

no  longer  lie  in  a  smooth  curve  as  originally  indicated  in 

2  4 

the  work  on  the  silicon  series.  Figure  4.12  summarises 
the  results  of  the  phosphorus  series,  the  transition 
potential  values  are  omitted  for  clarity.  All  the  new 
values  are  presented  in  Table  4.9. 

D.  Results  of  Calculations  -  Application  of  the  Various 

Models  and  Comparison  with  Experiment 

The  charges  and  potentials  obtained  by  methods  de¬ 
scribed  in  the  previous  section  can  be  used  in  the  "atom 
charge  model"  (Equation  3.1)  to  estimate  the  binding  energy 
shifts  of  the  core  electrons  emitted  upon  the  initial 
photo-ionisation  (AE(ls)  or  AE(2p)),  and  the  binding 

energy  shifts  of  a  2p  electron  in  a  system  already  con- 

* 

taining  a  2p  vacancy  (AE(2p  )).  The  latter  can  be  related 
to  the  Auger  energy  shifts  (AEtKL^L^))  by  subtraction  of 
the  Auger  shifts  from  the  x-ray  shifts  as  shown  in  Equation 

3.6. 

An  estimation  of  the  k  and  £  values  required  in 
Equation  3.1  may  be  obtained  from  a  least  squares  fit  of 


AE  -  V.=  kq  +  Z 


(4.2) 


Table  4.9.  Central  atom  changes  (q)  and  off-atom  potential  from  EWMO  calculation  (New  P  parameterisation). 
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State  of  Core  Ionisation 


Figure  4 . 12  As  Figure  4.6  but  with  new  P  parameterisation . 
The  transition  states  have  been  omitted  for  clarity. 
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where  AE  is  the  appropriate  experimental  chemical  shift, 
and  q  and  V  are  the  charges  and  potentials  appropriate  to 
the  model  under  consideration. 

The  ground  potential  model  (GPM)  estimates  AE(ls) 
and  AE(2p)  from  ground-state  charges  and  potentials 
(qN,  VN) .  The  charges  and  potentials  of  the  singly-ionised 
core-hole  state  (qN  ^ ,  VN  ^)  ,  obtained  by  using  the  (Z+l)+ 

equivalent  core  for  the  Z  atom,  provided  GPM  estimates  of 

* 

AE(2p  ).  The  relaxation  potential  model  (RPM)  estimates 

AE(ls)  and  AE(2p)  using  1/2  (qN+qN  and  1/2  (VN+VN  for  the 

*  N— 1  N— 2 

charges  and  the  potentials  and  AE(2p  )  from  1/2 (q  +q  ) 

and  1/2(V  +V  )  .  In  the  latter  case  q  and  V  are 

the  charges  and  potentials  of  the  doubly-ionised  core-hole 

states  obtained  by  using  the  (Z+2)  equivalent  core  for 

the  Z  atom. 

Both  the  GPM  and  RPM  approaches  have  been  used  in 

conjunction  with  charges  and  potentials  obtained  from 

CNDO/2  and  from  EWMO  (Hiickel)  calculations.  The  transition 

potential  model  (TPM) ,  which  uses  the  transition  charges 

and  potentials  (qN'N_1,  VN'N_1)  obtained  from  equivalent  core 

(Z+l/2 ) +1//2  calculations  was  applied  to  estimate  AE(ls) 

and  AE(2p).  AE(2p  )  was  estimated  using  (q  ,  V  ) 

+  3/2 

as  estimated  with  a  (Z+3/2)  7  equivalent  core  calculation. 

Only  charges  and  potentials  from  the  EWMO  (Hiickel)  calc¬ 
ulation  were  used  for  the  TPM  case. 
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A  theoretical  value  for  k  can  also  be  obtained  within 
the  context  of  the  point-charge  approximation 

=  < j| l/ri| =  n  (4-3) 

where  £  is  the  Slater  orbital  exponent  of  the  valence 
shell  and  n  is  the  principal  quantum  number  of  that  shell. 
The  value  for  £  is  unambiguous  within  the  CNDO/2  formalism, 
however,  within  the  EWMO  (Hiickel)  program  a  different  £  is 
used  for  the  s  and  p  valence  orbitals.  Table  4.10 
lists  the  appropriate  theoretical  k's  estimated  from  the 
Slater  exponents.  A  weighted  average  is  used  for  the 
EWMO  (Hiickel)  exponents.  The  Z+l/2  values  are  the  average 
of  the  exponents  appropriate  to  Z  (the  atomic  number)  and 
Z+l  atoms. 

The  GPM  approach  assumes  that  the  shifts  in  binding 

energy  can  be  connected  exclusively  to  the  properties  of 

the  ground  state  of  the  molecule.  In  effect  a  constant 

relaxation  contribution  within  a  series  of  molecules  is 

assumed.  Furthermore,  the  use  of  a  fitting  technique  to 

extract  predictive  parameters  (k  and  Z)  from  experimental 

data  tends  to  obscure  differences  arising  from  additional 

contributions  which  may  be  significant.  Incorporation  of 

relaxation,  as  in  the  use  of  the  relaxation  or  transition 

potential  models,  should  give  a  more  correct  estimation  of 

* 

k.  The  values  for  k  obtained  by  fitting  the  2p  and  2p 
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Table  4.10.  Theoretical  k' 

1 sa  from  Slater's 

exponents . 

Equivalent 

P 

Core 

S 

EWMOb 

CNDO/2 

P 

15.3446 

14.5124 

16.2938 

15.4950 

s+ 

S 

17.2430 

16.4776 

18.2271 

17.4601 

Cl++ 

Cl+ 

19.2112 

18.4425 

20.2061 

19.4257 

Ar++ 

21.2010 

20.4080 

(a)  k  =  C/3  x  27.2107  ev. 
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experimental  shifts  with  RPM  (or  TPM)  charges  and  potentials 
should  display  a  similar  trend  to  the  theoretical  values. 
This  relationship  provides  an  indication  of  whether  the 
charges  and  potentials  obtained  from  the  calculations  truly 
reflect  the  differences  between  the  photo-ionisation  of 
the  molecule  and  the  loss  of  an  electron  from  the  singly- 
charged  ion . 

(i)  CNDO/2  Calculations 

The  values  obtained  for  k  and  Z  from  the  least-squares 
fit  are  listed  in  Table  4.11.  The  GPM  and  RPM  shifts  for 
the  phosphorus  series  are  shown  in  Table  4.12  and  those 
for  sulphur  in  Table  4.13.  The  correlations  obtained  be¬ 
tween  the  various  calculated  shifts  and  experiment  are 
summarised  in  Table  4.14.  The  correlations  obtained  between 
the  RPM  calculated  shifts  and  experiment  for  the  2p  and  Z 
Auger  shifts  are  shown  in  Figures  4.13  to  4.16. 

In  general  the  models  provide  reasonable  correlation 
with  experiment.  The  RPM  approach  gives  a  slope  with  almost 
the  ideal  value  of  unity  in  all  cases  (a  suitable  choice  of 
reference  would  give  zero  intercept) .  The  phosphorus  series 
shows  clearly  the  improvement  in  correlation  resulting  from 
the  RPM  approach  as  compared  to  the  GPM  approach.  The  im¬ 
provement  is  less  pronounced  for  the  sulphur  series.  The 
correlations  obtained  for  the  Auger  data  for  both  the 
phosphorus  and  sulphur  series  were  not  as  good  as  those 
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Table  4.11.  Parameters  from  least-squares  fita  for  the  atom- 

charge  model  (CNDO/2  charges  and  potentials) . 


k 

£ 

corr . 
coef f . 

Phosphorus  series 

Is  (GPM) 

14.2372 

-  1.9357 

0.9949 

Is  (RPM) 

14.2604 

-  5.6186 

0.9988 

2p  (GPM) 

13.6076 

-  1.9476 

0.9944 

2p  (RPM) 

13.5715 

-  5.5893 

0.9988 

2p*  (GPM) 

14.6411 

-12.1232 

0.9940 

2p*  (RPM) 

14.3854 

-15.3313 

0.9914 

Sulphur  series 

Is  (GPM) 

15.3715 

-  0.2772 

0.9980 

Is  (RPM) 

15.8820 

-  4.8491 

0.9977 

2p  (GPM) 

14.6091 

-  0.4162 

0.9980 

2p  (RPM) 

14.9901 

-  4.8689 

0.9978 

2p*  (GPM) 

15.5111 

-11.6618 

0.9964 

2p*  (RPM) 

16.1632 

-16.2604 

0.9890 

(a)  Data  fitted  to 

AE-V  =  kq  +  £ 
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Table  4.12.  Calculated  atom  charge  model  shifts  (ev)  obtained  from  CNDO/2  charges  and 
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Correlation  obtained  from  least-squares  fit  on 
calculated  versus  experimental  shifts  (CNDO/2 
calculation : atom-charge  model). 

corr . 


slope 

intcpt 

coef f . 

Phosphorus  series 

Is  (GPM) 

0.9522 

1.3943 

0.9504 

Is  (RPM) 

1.0006 

-0.6136 

0.9903 

2p  (GPM) 

0.9377 

1.5246 

0.9353 

2p  (RPM) 

1.0003 

-0.4484 

0.9884 

KL2L3  (GPM) 

0.9466 

0.2283 

0.9615 

KL2L3  (RPM) 

1.0346 

-1.6707 

0.9629 

Sulphur  series 

Is  (GPM) 

0.9908 

0.6749 

0.9847 

Is  (RPM) 

0.9972 

-1.0617 

0.9864 

2p  (GPM) 

0.9911 

0.7469 

0.9818 

2p  (RPM) 

0.9999 

-0.9586 

0.9848 

KL2L3  (GPM) 

0.9997 

-0.4783 

0.9765 

KL2L3  (RPM) 

1.0250 

-2.1795 

0.9555 

100 


AE0  (expt)/cV 
2p 


Figure  4.13  Correlation  of  calculated  versus  experimental 
P2p  binding  energy  shifts.  Calculated  shifts  are  ob¬ 
tained  from  the  atom  charge  model  using  CNDO/2  RPM  charges 


and  potentials. 
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Figure  4.14  Correlation  of  calculated  versus  experimental 
S2p  binding  energy  shifts.  Calculated  shifts  are  obtained 
from  the  atom  charge  model  using  CNDO/2  RPM  charges  and 
potentials . 
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Figure  4.15  Correlation  of  calculated  versus  experimental 
P  KL2L3  Auger  energy  shifts.  Calculated  shifts  are  ob¬ 
tained  from  the  atom  charge  model  using  CNDO/2  RPM  charges 


and  potentials. 
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Figure  4.16  Correlation  of  calculated  versus  experimental 
S  KI^L^  Auger  energy  shifts.  Calculated  shifts  are  ob¬ 
tained  from  the  atom  charge  model  using  CNDO/2  RPM  charges 
and  potentials. 
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obtained  for  the  core-electron  shifts.  This  suggests 
limitations  in  the  applicability  of  the  equivalent  core 
approximation  for  the  doubly-charged  ion,  especially  with 
charges  and  potentials  obtained  by  the  CNDO/2  methodology. 

The  value  obtained  for  k  for  the  2p  (RPM)  shift  is  larger 
than  that  for  the  2p  (RPM)  shift,  however,  neither  value 
is  as  large  as  the  theoretical  values.  This  suggests  that 
the  calculation  is  able  to  predict  the  general  trend  in  a 
satisfactory  manner  although  the  charges  and  potentials 
may  not  be  absolutely  correct. 

The  shifts  obtained  by  the  potential  model  (equation  3.3) 
have  also  been  evaluated  and  the  results  are  presented  in 
Table  4.15.  The  correlations  between  the  calculated  and 
experimental  shifts  are  summarised  in  Table  4.16  and  those 
from  the  relaxation  potential  model  are  shown  in  Figures 
4.17  to  4.22.  The  correlations  obtained  between  calculated 
and  experimental  shifts  are  as  good  or  better  than  those 
obtained  with  the  atom-charge  model,  however,  the  slope 
of  the  line  is  no  longer  unity.  In  this  model  a  theor¬ 
etical  value  for  k  is  calculated  each  time  using  k^ = < j | 1/r^ | j > 
(equation  3.2).  This  supports  the  point  that  the  valence 
populations  (and  hence  charges)  and  potentials  obtained  in 
the  CNDO/2  calculations  are  not  consistent  with  the 
theoretically  estimated  values  for  k  but  are  still  adequate 
in  showing  the  trends. 
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TABLE  4.15  Calculated  potential  model  shifts  (eV)  obtained  from  CNDO/2  calculation 
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Table  4.16.  Correlation 

obtained 

from  least-squares  fit  on 

calculated 

versus  experimental  shifts 

(CNDO/2 

calculation 

rpotential 

model) . 

corr . 

slope 

intcpt 

coef f . 

Phosphorus  series 

Is  (GPM) 

1.0066 

1.4230 

0.9527 

Is  (RPM) 

1.2124 

-0.6626 

0.9914 

2p  (GPM) 

1.1371 

1.6727 

0.9427 

2p  (RPM) 

1.3767 

-0.3814 

0.9862 

KL2L3  (GPM) 

1.2158 

1.0925 

0.9593 

KL2L3  (RPM) 

1.4372 

-0.5510 

0.9761 

Sulphur  series 

Is  (GPM) 

1.1869 

0.7235 

0.9878 

Is  (RPM) 

1.2122 

-1.1111 

0.9895 

2p  (GPM) 

1.3736 

0.8730 

0.9876 

2p  (RPM) 

1.4049 

-0.9644 

0.9908 

KL2L3  (GPM) 

1.4662 

0.6209 

0.9845 

KL2L3  (RPM) 

1.4080 

-1.3835 

0.9766 
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Figure  4.17  Correlation  of  calculated  (potential  model) 
shifts  versus  experimental  Pis  binding  energy  shifts.  Cal¬ 
culated  shifts  are  obtained  using  CNDO/2  RPM  point  charge 
potentials . 
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Figure  4.18  Correlation  of  calculated  (potential  model) 
shifts  versus  experimental  Sis  binding  energy  shifts. 
Calculated  shifts  are  obtained  using  CNDO/2  RPM  point 


charge  potentials. 
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Figure  4.19  Correlation  of  calculated  (potential  model) 
shifts  versus  experimental  P2p  binding  energy  shifts. 
Calculated  shifts  are  obtained  using  CNDO/2  RPM  point 
charge  potentials. 
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Figure  4.20  Correlation  of  calculated  (potential  model) 
shifts  versus  experimental  S  binding  energy  shifts. 
Calculated  shifts  are  obtained  using  CNDO/2  RPM  point 
charge  potentials. 
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Figure  4.21  Correlation  of  calculated  (potential  model) 
shifts  versus  experimental  P  KL^L^  Auger  energy  shifts. 
Calculated  shifts  are  obtained  using  CNDO/2  RPM  point 
charge  potentials. 
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Figure  4.22  Correlation  of  calculated  (potential  model) 
shifts  versus  experimental  S  KL^L^  Auger  energy  shifts. 
Calculated  shifts  are  obtained  using  CNDO/2  RPM  point 
charge  potentials. 
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The  slope  obtained  for  the  Is  shift  prediction  is 
closer  to  unity  than  that  for  the  equivalent  2p  shift, 
reflecting  the  fact  that  the  point-charge  approximation 
is  more  appropriate  for  the  deeper  levels.  The  2p  electrons 
will  have  a  greater  interaction  with  the  valence  shell  than 
will  the  Is  level  because  of  its  closer  proximity. 

CNDO/2  calculations  with  d-orbitals  included  in  the 
basis  set  were  performed  for  the  phosphorus  series.  Con¬ 
vergence  problems  occurred  in  many  cases  and  so  only  the 
results  of  a  limited  series  are  presented.  Table  4.17 
shows  the  calculated  shifts  along  with  the  correlations 
obtained  from  least-squares  fits  on  the  calculated  versus 
experimental  shifts.  The  correlations  obtained  are  con¬ 
siderably  poorer  than  those  obtained  from  calculations 
omitting  d-orbitals.  Furthermore,  the  values  for  k  (Table 
4.18),  though  larger  than  those  obtained  when  d-orbitals 
were  omitted,  do  not  follow  the  expected  theoretical  trend 
when  comparing  the  loss  of  a  2p  electron  from  the  molecule 
to  the  loss  of  a  2p  electron  from  the  ion.  It  can  be 
concluded  that  the  charges  and  potentials  obtained  from 
the  CNDO/2  calculation  including  d-orbitals  in  the  basis 
set  are  inadequate. 

(ii)  EWMO  (Huckel )  Calculations 

Only  the  atom-charge  model  was  used  in  conjunction  with 
charges  and  potentials  from  the  EWMO  (Huckel)  calculation. 


<n 


j  I 


TABLE  4.17  Calculated  atom  charge  model  shifts  (eV)  obtained  from  CNDO/2  charges 
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Table  4.18.  k  values  obtained  from  least-squares  fita  using 

CNDO/2  charges  and  potentials  -  d  orbitals  in¬ 
cluded  in  the  basis  set. 


Phosphorus  series  k 

Is  (GPM)  18.9895 

Is  (RPM)  17.6137 

2p  (GPM)  17.8466 

2p  (RPM)  16.4861 

2p*  (GPM)  17.8507 

2p*  (RPM)  16.1429 


(a)  k  obtained  from  fit  to 

AE-V  =  kq  +  l 


. 
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The  fitted  values  of  k  and  i  for  the  GPM,  RPM  and  TPM 
approaches  are  listed  in  Table  4.19.  The  results  of 
calculations  done  with  the  original  parameterisation  and 
with  new  phosphorus  parameters  (Appendix  2)  are  compared. 

The  shifts  for  the  phosphorus  series  are  given  in  Table 
4.20  for  the  case  of  the  original  parameterisation  for 
P  and  in  Table  4.21  for  the  new  parameters.  Shifts  for 
the  sulphur  series  are  given  in  Tables  4.22  and  4.23 
respectively.  Table  4.24  summarises  the  correlations  be¬ 
tween  the  various  calculated  shifts  and  experiment,  some 
of  which  are  shown  in  Figures  4.23  to  4.30. 

The  correlations  between  the  calculated  and  experimental 
shifts  are  very  similar  for  both  parameter isations .  The 
original  parameterisation  gives  a  slightly  better  fit  for 
the  phosphorus  series  whereas  the  converse  is  true  for  the 
sulphur  series  in  which  a  phosphorus  group  is  present  as 
a  substituent  on  the  sulphur.  The  use  of  fitting  tech¬ 
niques  to  obtain  k  and  i  values  tends  to  obscure  the 
effect  that  the  new  parameterisation  may  have  on  the  fit. 

The  new  parameterisation  clearly  gives  a  lower  value  for 
the  k  associated  with  the  initial  photo-ionisation  in  the 
phosphorus  series  and  this  may  indicate  that  the  new 
parameterisation  is  more  realistic,  however  it  is  difficult 
to  draw  any  definite  conclusions. 

The  correlations  between  experimental  and  calculated 
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Table  4.19.  Parameters  from  least-squares  fit  for  atom-charge  model  (EWMO  charges 
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Table  4.20.  Calculated  atom  charge  model  shifts  (ev)  obtained  from  EWMO  (original  parameterisation)  charges  and 
potentials3'^  -  Phosphorus  series. 
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Table  4.21.  Calculated  atom  charge  model  shifts  (eV)  ob¬ 
tained  from  EWMO  (new  P  parameterisation) 
charges  and  potentials  useda. 


L2p 


MOLECULE 

GPM 

RPM 

ph3 

0.0 

0.0 

spf3 

7.18 

6.10 

p(ch3)3 

0.48 

-0.82 

spci3 

6.26 

4.73 

opf3 

7.70 

6.85 

0PC13 

6.08 

4.59 

PFC 

5 

9.12 

8.09 

PCI3 

3.05 

1.61 

PF3 

5.44 

4.58 

SP  (CH30)  Cl 

4.90 

3.41 

(CH30) 3P0 

4.20 

2.89 

(CH30) 3PS 

4.11 

2.71 

(ch30)3p 

1.50 

0.24 

P  (CH3)C12 

2.29 

0.89 

OP  (CH2C1)C12 

5.77 

4.16 

SP (CH3)C12 

5.68 

4.15 

P(C2H5}3 

0.52 

-1.24 

(a)  2p*  and  KL^L^  shifts  are 


TPM 

GPM 

RPM 

TPM 

0.0 

0.0 

0.0 

0.0 

5.92 

6.38 

5.33 

5.16 

-1.05 

0.43 

-0.84 

-1.06 

4.51 

5.77 

4.28 

4.07 

6.69 

6.73 

5.91 

5.75 

4.39 

5.43 

3.98 

3.78 

7.91 

8.07 

7.07 

6.89 

1.37 

2.84 

1.45 

1.21 

4.42 

4.83 

4.01 

3.85 

3.19 

4.24 

2.80 

2.58 

2.67 

3.30 

2.03 

1.81 

2.48 

3.36 

2.00 

1.78 

-0.01 

0.99 

-0.22 

-0.46 

0.66 

2.12 

0.77 

0.54 

3.95 

5.14 

3.57 

3.37 

3.93 

5.21 

3.72 

3.50 

-1.49 

0.49 

-1.21 

-1.46 
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Table  4.23.  calculated  atom  charge  model  shifts  (eV)  obtained  from  EWMO  (new  P  parameterisation)  charges 
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(a)  The  shifts  obtained  for  the  non-P  containing  compounds  may  be  different  from  those  given  in  Table  4.22  since 


TABLE  4.24  Correlation  obtained  from  least-squares  fit  on 


calculated  versus  experimental  shifts  (EWMO 
calculation) 

Original  Parametrisation  New  P  Parametrisation 


slope 

intcpt 

corr  coeff 

slope 

intcpt 

corr  coeff 

Phosphorus  Series 

Is  (GPM) 

0.9649 

0.7501 

0.9617 

0.9626 

1.0972 

0.9490 

Is  (RPM) 

0.9802 

-0.4154 

0.9682 

0.9800 

-0.2174 

0.9658 

Is  (TPM) 

0.9819 

-0.5537 

0.9672 
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Figure  4.23  Correlation  of  calculated  versus  experimental 
P2p  binding  energy  shifts.  Calculated  shifts  are  obtained 
from  the  atom  charge  model  using  EWMO  RPM  charges  and 


potentials . 
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AE2p (expt) /eV 


Figure  4.24  Correlation  of  calculated  versus  experimental 
S2p  binding  energy  shifts.  Calculated  shifts  are  obtained 
from  the  atom  charge  model  using  EWMO  RPM  charges  and 
potentials . 
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Figure  4.25  Correlation  of  calculated  versus  experimental 


P2p  binding  energy  shifts.  Calculated  shifts  are  obtained 
from  the  atom  charge  model  using  EWMO  (new  P  parameterisation) 
RPM  charges  and  potentials. 


. 


130 


-50  5  10 

AE2p  (cxpt)  /eV 


Figure  4.26  Correlation  of  calculated  versus  experimental 
S2p  binding  energy  shifts.  Calculated  shifts  are  obtained 
from  the  atom  charge  model  using  EWMO  (new  P  parameterisation ) 
RPM  charges  and  potentials. 
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Figure  4.27  Correlation  of  calculated  versus  experimental 
P2p  binding  energy  shifts.  Calculated  shifts  are  ob¬ 
tained  from  the  atom  charge  model  using  EWMO  (new  P 
parameterisation)  TPM  charges  and  potentials. 
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Figure  4.28  Correlation  of  calculated  versus  experimental 


S2p  binding  energy  shifts.  Calculated  shifts  are  obtained 
from  the  atom  charge  model  using  EWMO  (new  P  parameter- 
isation)  TPM  charges  and  potentials. 


. 


1  0 


5 


> 

o 

1  0 

H 


-5 


-  1  0 

-  1  0  -5  0  5  1 

AEKLL (exPt)/eV 


Figure  4.29  Correlation  of  calculated  versus  experimental 
P  k^2L3  Au9er  energy  shifts.  Calculated  shifts  are  ob¬ 
tained  from  the  atom  charge  model  using  EWMO  RPM  charges 
and  potentials. 
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Figure  4.30  Correlation  of  calculated  versus  experimental 
S  KL^L^  Auger  energy  shifts.  Calculated  shifts  are  ob¬ 
tained  from  the  atom  charge  model  using  EWMO  (new  P  para- 
meterisation)  RPM  charges  and  potentials. 
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shifts -for  both  the  RPM  and  TPM  approaches  are  equally 
good.  The  TPM  approach  gives  a  slightly  lower  value  for  k. 
Since  the  TPM  approach  was  not  used  with  the  CNDO/2  formal¬ 
isation  we  cannot  compare  the  advantages  of  the  TPM  and 
thus  we  shall  only  consider  the  RPM  approach.  However,  it 
seems  that  the  TPM  approach  gives  as  reasonable  an  estimate 
of  chemical  shifts  at  this  level  of  approximation  as  does 
the  RPM  approach  and  does  so  with  one  calculation  rather 
than  two . 

The  EWMO  (Huckel)  calculations,  for  the  sulphur  series, 

gave  as  good  a  correlation  of  experimental  versus  calculated 

shifts  at  all  levels  as  the  CNDO/2  calculations,  however 

* 

the  value  for  k  obtained  on  the  RPM  fit  on  the  2p  data 

was  virtually  the  same  as  that  for  the  2p  data.  In  the 

phosphorus  series  the  correlations  derived  from  the  EWMO 

(Huckel)  calculations  were  considerably  worse  than  those 

derived  from  the  CNDO/2  calculations.  While  the  experimental 

versus  calculated  shift  correlations  for  the  Is  and  2p 

electrons  were  reasonable,  the  correlation  of  the  Auger 

data  was  poor.  The  model,  utilising  EWMO  (Huckel)  charges 

and  potentials,  exhibits  deficiencies  when  dealing  with  the 

loss  of  the  second  electron  and  the  concomitant  electron 

reorganisation.  This  is  demonstrated  by  the  significantly 

lower  correlation  coefficient  obtained  on  the  fit  of 

* 

equation  4.2  for  the  P2p  shift.  The  poor  fit  here  will 
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result  in  a  lack  of  correlation  between  the  experimental 
and  calculated  shifts.  Figure  4.31  shows  the  correlation 

"fc 

of  (AE  -  V)  with  q  for  the  P2p  shifts  using  the  RPM 

charges  and  potentials  appropriate  to  the  ion.  It  is  clear 

why  a  poor  correlation  was  obtained.  The  EWMO  (Huckel) 

calculation  results  for  phosphorus  compounds  with  a  lone 

pair  on  the  central  atom  are  different  from  those  without 

a  lone  pair.  This  difference  is  also  illustrated  by  a 

* 

similar  plot  for  the  S2p  shifts  (Figure  4.32).  The 

molecules  can  be  grouped  according  to  the  number  of  lone 

pairs  on  the  central  atom.  The  "pivoting"  is  not  as 

dramatic  in  the  sulphur  case  as  with  the  phosphorus  series 

and  hence  a  fit  over  all  the  groups  gives  a  reasonable 

correlation,  however,  the  pivoting  results  in  an  unreasonable 

value  for  k.  Similar  differences  in  the  correlation  of 

(AE  -  V)  against  q  are  observed  for  the  Pis  shifts  using 

the  RPM  charges  and  potentials  appropriate  to  the  molecule 

(Figure  4.33).  The  differences  are  not  as  pronounced  as 

* 

those  associated  with  the  P2p  shifts,  and  as  in  the 
sulphur  case,  results  in  a  reasonable  fit.  Correlation  of 

•k 

(AE  -  V)  with  q  for  the  P2p  shifts,  using  RPM  charges 
and  potentials  appropriate  to  the  ion  from  the  CNDO/2 
calculation  (Figure  4.34)  shows  a  consistent  trend  re¬ 
gardless  of  coordination. 

Table  4.25  summarises  the  k  and  £  values  obtained  with 
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(new  P  parameterisation)  RPM  charges  and  potentials  ap¬ 


propriate  to  the  ion  state  -  phosphorus  series.  Compound 
numbers  are  presented  in  Table  4.1.  The  compounds  are 
grouped  according  to  those  with  one  lone  pair  (Ki)  and  no 
lone  pairs  (©)  on  the  phosphorus  atom. 
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Figure  4.32  Correlation  of  (AE2p*-V)  with  (kq+£)  using  EWMO 
(new  P  parameterisation)  RPM  charges  and  potentials  appro¬ 
priate  to  the  ion  state  -  sulphur  series.  Compound  numbers 
are  presented  in  Table  4.2.  The  compounds  are  grouped  ac¬ 
cording  to  those  with  three  lone  pairs  (V)  ,  two  lone  pairs  (□)  , 
one  lone  pair  (A)  and  no  lone  pairs  (0)  on  the  sulphur  atom. 
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Figure  4.33  Correlation  of  (AE1 „-V)  with  (kq+£)  using  EWMO 

"  _L  S 

(new  P  parameterisation)  RPM  charges  and  potentials  ap¬ 
propriate  to  the  molecule  -  phosphorus  series.  Compound 
numbers  are  presented  in  Table  4.1.  The  compounds  are 
grouped  according  to  those  with  one  lone  pair  (O)  and  no 
lone  pairs  (□)  on  the  phosphorus  atom. 
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Figure  4.34  Correlation  of  (AE2p*~V)  with  (kq+£)  using 
CNDO/2  RPM  charges  and  potentials  appropriate  to  the  ion 
state  -  phosphorus  series.  Compound  numbers  are  presented 


in  Table  4.1. 
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EVJMO  (Hiickel)  charges  and  potentials,  as  well  as  the  cor¬ 
relation  between  the  experimental  and  calculated  shifts, 
for  the  phosphorus  (original  and  new  parameterisation )  and 
sulphur  series  in  which  only  molecules  without  lone  pairs 
on  the  central  atom  have  been  considered.  The  correlation 
obtained  is  a  clear  improvement  over  that  when  all  the 
molecules  are  considered.  The  k  values  obtained  for  the 
RPM  approach  now  reflect  the  theoretical  trend.  The  new 
parameterisation  gives  lower  values  for  k  but  again  the 
differences  are  too  small  to  permit  definite  conclusions. 

The  EWMO  (Hiickel)  calculation  seems  to  partition 
charge  differently  when  lone  pairs  are  involved  on  the  atom 
and  as  a  result  series  of  molecules  related  by  equal 
numbers  of  lone  pairs  are  established.  It  must  be  concluded 
that  the  general  applicability  of  the  EWMO  (Hiickel) 
calculation  to  estimate  shifts  is  limited  except  in  the 
case  of  a  carefully  selected  series  of  compounds. 

E.  Relationship  of  the  Auger  Electron  Shift  to  the  Core 

Electron  Binding  Energy  Shift 

34 

The  Auger  parameter  was  originally  defined  as 

Aa  =  (AEauq  +  AEbe)  =  2ARea  (4.4) 

The  sum  of  the  Auger  chemical  shift  and  the  binding  energy 
shift  of  a  specified  core  electron  equals  twice  the  extra 
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atomic  -relaxation  energy  associated  with  the  photoelectron. 
It  was  shown  in  the  introduction  that  a  more  accurate 
expression  can  be  obtained.  For  the  KLL  Auger  shifts  the 
Auger  parameter  becomes 


L+K 

X-ray 


(AE„tt  +  AE  )  =  ARt  +  AR*  +  AE 

1a±j1j  J_i  J_j  J_j 


(4.5) 


with  the  specified  core  electron  coming  from  the  L-shell. 

The  Auger  parameter  depends  not  only  on  the  relaxation 
energy  associated  with  the  loss  of  the  initial  photoelectron 
but  also  with  that  of  the  second  electron  which  is  leaving 
a  singly-charged  ion.  The  X-ray  term  corrects  for  the 
difference  in  the  chemical  shift  associated  with  the  K-shell 
and  that  associated  with  the  L-shell,  this  must  be  considered 
since  the  K-shell  is  refilled  by  an  L-electron  in  the  Auger 
process.  If  the  specified  electron  is  from  the  K  shell, 
the  Auger  parameter  becomes 


(4.6) 


(iEKLL  +  AV  =  ARL  +  ARL*  +  2AEX-ray 


The  relaxation  terms  are  still  associated  with  the  L  shell, 
the  shell  which  provides  the  Auger  electron  which  we  observe. 

The  values  for  AR  in  the  above  equations  express 
changes  in  total  relaxation  energies.  If  the  atomic 
contribution  is  considered  to  be  constant  for  the  whole 
series  then  the  values  for  AR  are  shifts  in  extra  atomic 


relaxation . 
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' 
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Rearrangement  of  equation  (4.5)  gives 


AEKLL  ^ARL  +  ARL*}  “  ^EL  “  AEX-ray ^  (4.7) 

The  Auger  chemical  shift  depends  on  a  reference  level 
ir^K 

(AEl  -  ^Ex-ray^  w^ich  is  modified  by  the  relaxation  con¬ 
tributions  arising  from  the  removal  of  both  the  initial  and 
secondary  electrons.  This  expression  provides  a  convenient 
way  of  relating  the  Auger  shift  to  the  binding-energy  shift. 
This  is  shown  diagramatically  in  Figure  4.35  for  the 
phosphorus  series  and  Figure  4.36  for  the  sulphur  series. 

Only  experimental  values  have  been  used  and  hence  the 
relaxation  contribution  shown  is  the  sum  of  AR  and  AR  * 

1j  Xj 

contributions.  The  diagram  clearly  indicates  where  and  why 
reversals  occur  between  the  Auger  and  photoelectron  chemical 
shifts.  The  relaxation  contribution  can  greatly  alter  the 
shifts — adding  as  much  as  3.5  eV  in  some  cases  relative 
to  the  reference. 

The  results  from  the  calculations  have  also  been 
applied  to  equation  4.5.  In  view  of  the  discussion  in  the 
previous  section  only  the  results  from  the  CNDO/2  calculation 
will  be  considered  for  the  complete  series  of  molecules. 

The  Auger  parameter,  experimental  relaxation  energy  shifts, 
as  well  as  the  calculated  relaxation  energy  shifts  for  both 
electrons  obtained  from  the  atom  charge  and  potential 
models  are  shown  in  Table  4.26  for  the  phosphorus  series 
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Figure  4.35  Relationship  of  the  Auger  energy  shift  to  the 
photoelectron  binding  energy  shift  for  the  phosphorus 
series.  Compound  numbers  are  presented  in  Table  4.1. 
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Figure  4.36 


Relationship  of  the  Auger  energy  shift  to  the 


photoelectron  binding  energy  shift  for  the  sulphur  series. 
Compound  numbers  are  presented  in  Table  4.2. 


tabu;  4.26  Relaxation  -  experimental  and  from  CNDO/2 

calculation  and  the  Auger  parameter  -  P  series 
(all  shifts  in  eV) 


Atom  Charge  Modelc  Potential  Model 


Molecule 

Total 

Relax 

aexray 

— - - . 

Aa  b 

arl 

arl+arl 

+aexray 

arl 

arl 

AVARi 

+aexray 

ph3 

0.00 

0 

0 

0 

0 

0 

0 

0 

0 

spf3 

1.12 

1.05 

2.17 

0.67 

0.86 

2.58 

0.33 

0.76 

2.14 

P(CH3}3 

2.08 

0.12 

2.20 

1.68 

1.61 

3.41 

1.77 

1.56 

3.45 

spci3 

3.40 

0.64 

4.04 

1.95 

1.99 

4.58 

1.99 

2.16 

4.79 

opf3 

-0.28 

1.05 

0.77 

0.15 

0.34 

1.54 

-0.33 

0.10 

0.82 

0PC13 

2.60 

0.69 

3.29 

1.73 

1.95 

4.37 

1.68 

2.10 

4.47 

pf5 

0.18 

1.22 

1.40 

0.52 

0.78 

2.52 

0.07 

0.72 

2.01 

PC13 

2.69 

0.49 

3.18 

1.68 

1.98 

4.15 

1.70 

2.17 

4.36 

PF2 

-0.66 

0.76 

0.10 

0.12 

0.39 

1.27 

-0.34 

0.11 

0.53 

SP(Ch30)2Cl 

2.91 

0.82 

3.73 

2.10 

2.09 

5.01 

2.13 

2.26 

5.21 

<ch30) 3po 

1.85 

0.80 

2.65 

1.89 

1.97 

4.66 

1.82 

2.09 

4.71 

(Cfa'30)  3PS 

2.47 

0.87 

3.34 

2.11 

2.04 

5.02 

2.11 

2.18 

5.16 

(ch30) 3P 

1.59 

0.63 

2.22 

1.78 

1.89 

4.30 

1.75 

2.00 

4.38 

P(CH3)C12 

2.50 

0.36 

2.86 

1.65 

1.84 

3.85 

1.67 

1.94 

3.97 

p (C2^5 ) 3 

2.90 

0.05 

2.95 

2.29 

2.19 

4.53 

2.53 

2.32 

4.90 

OP(Ch_Cl)Cl2 

2.68 

0.59 

3.28 

1.93 

1.97 

4.49 

1.92 

2.06 

4.57 

SP(CH3)C12 

3.43 

0.60 

4.03 

1.96 

1.91 

4.47 

2.01 

2.03 

4.64 

P(CF3)3 

2.82 

0.04 

2.86 

2.19 

2.06 

4.29 

2.48 

2.21 

4.73 

(a)  (AR^ARjO  =  Aa-AEXRAY 
(h)  Aa  =  AEKLli  +  AEL 

(c)  k's  used  GPM(2p)  =  13.1646  RPM(2p)  =  13.5715 

GPM (2p*)  =  13.9785  RPM(2P*)  =  14.3854 

(d)  Correlation  with  Aa  (see  Fig.  4.37) 

slope  =  1.0899  inept  =  0.8845  corr  coeff  =  0.9281 

(e)  Correlation  with  Aa  (see  Fig.  4.37) 
slope  =  1.2969  inept  =  0.3549  corr  coeff  = 


0.9323 
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and  in  Table  4.27  for  the  sulphur  series. 

The  ground  potential  values  for  k  used  in  the  atom- 
charge  model  were  obtained  from  a  linear  extrapolation 
based  upon  the  relaxation  potential  values  for  k.  Figures 
4.37  and  4.38  show  the  correlation  of  (AE  +  AE  )  with 

KL»Ju  -  Xj 

Lr>- 

(AR^  +  ARl*  +  AEx-ray^  f°r  Phosphorus  and  sulphur 

series  respectively.  The  correlation  is  satisfactory  with 
the  sulphur  series  having  less  scatter  than  the  phosphorus 
series.  In  both  cases  the  potential  model  (utilising 
theoretical  k's)  gives  the  better  correlation.  Unit  slopes 
are  not  obtained,  which  is  not  surprising  in  view  of  the 
correlations  obtained  between  the  experimental  and  cal¬ 
culated  shifts  presented  previously.  With  the  exception 
of  the  molecules  containing  the  (-OCH^)  ligand,  the  cal¬ 
culated  values  follow  the  same  general  trends  as  experiment 
though  they  are  somewhat  larger. 

The  CNDO/2  calculation  used  here  fails  with  (-OCH3) . 

54 

Lee  et  al.  have  calculated  the  "electron  flow"  (relaxation) 
during  core-ionisation  for  a  series  of  tervalent  phosphorus 
compounds  using  a  CNDO/2  program  which  was  modified  by 
Sherwood.55  They  found  that  upon  core  ionisation,  the 
P (OCH^ ) 3  molecule  exhibited  substantially  less  relaxation 
than  PCI 3  and  P(CH3)3. 

The  results  from  the  EWMO  (Huckel)  calculations  have 
been  applied  to  select  series  of  compounds.  Each  series 
contains  molecules  with  the  same  number  of  lone  pairs  on  the 
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central  atom.  The  theoretical  k's  (Table  4.10)  were  used 
in  conjunction  with  the  atom-charge  model.  The  results, 
along  with  the  CNDO/2  relaxation  energy  shifts  based  upon 
the  same  reference  values,  are  shown  in  Table  4.28  for 
phosphorus  and  Table  4.29  for  sulphur.  The  trends  in  the 
total  relaxation  obtained  from  the  EWMO  (Hiickel)  cal¬ 
culations  follow  the  experimental  trends  in  all  cases 
(except  PF^  compared  to  PH^) .  The  values  obtained  are 
smaller  than  those  obtained  by  the  CNDO/2  method  and  the 
total  relaxation  energy  shifts  are  in  closer  agreement  with 
experiment.  In  some  cases  the  relative  relaxation  con¬ 
tributions  of  the  first  and  second  ionisation  steps  to  the 
total  relaxation  have  been  reversed  over  those  predicted 
by  the  CNDO/2  calculations. 

In  view  of  the  difficulties  and  contradictions  arising 
from  the  semi-empirical  calculations  used  here,  it  is 
difficult  to  draw  any  definite  conclusions  about  the  sig¬ 
nificance  of  changes  in  relaxation  energies  which  are 
associated  with  the  first  ionisation  step  compared  to  those 
associated  with  the  second.  However,  certain  observations 
can  be  made  based  upon  the  "experimental"  values  for  total 
relaxation.  The  relaxation  depends  on  the  ability  of  a 
group  to  release  electrons.  This  may  be  due  to  its  electro¬ 
negativity,  polarisability  or  both.  With  all  things  being 
equal  the  more  electronegative  a  substituent  the  less 
relaxation  will  occur.  This  is  amply  illustrated  by  the 
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Table  4.28.  Experimental  and  calculated  (atom  charge  model) 


relaxation 

energy  shifts 

(eV)  -  P 

.  b 
series 

• 

Exptl . 

CNDO/2 

EWMO 

Molecule 

Ar  +Ar 

Xj  Ij 

,*  iRL 

Arl* 

Arl 

Arl* 

PK3 

0 

0 

0 

0 

0 

PF3 

-0.66 

-0.34 

0.11 

0.14 

0.35 

PC13 

2.69 

1.70 

2.17 

1.33 

1.65 

P  (CH3)  3 

2.08 

1.77 

1.56 

1.35 

1.36 

P (OCH3) 3 

1.59 

1.75 

2.00 

0.75 

0.81 

P(CH3)C12 

2.50 

1.67 

1.94 

1.32 

1.56 

P(C2h5)3 

2.90 

2.53 

2.32 

1.92 

— 

P  (CF3)  3 

2.82 

2.48 

2.21 

— 

— 

PFC 

D 

0 

0 

0 

0 

0 

SPF3 

0.94 

0.26 

0.04 

0.42 

0.43 

SPC13 

3.22 

1.96 

1.44 

1.33 

1.39 

SP (OCH3) 

2.29 

2.04 

1.46 

0.83 

0.77 

SP  (OCH3)  Cl 

2.73 

2.06 

1.54 

1.05 

1.07 

SP  (CH3)C12 

3.45 

1.94 

1.31 

1.34 

1.25 

OPE  3 

-0.46 

-0.40 

-0.62 

00 

0 

• 

0 

1 

-0.15 

0PC13 

2.42 

1.61 

1.38 

1.06 

1.10 

OP  (CH30) 3 

1.67 

1.75 

1.37 

0.53 

0.44 

OP  (CH2C1)C12 

2.50 

1.85 

1.34 

1.22 

1.26 

(a)  Theoretical 

values 

for  k  used  (Table  4.10) 

• 

(b)  The  molecules  have 

been  grouped 

according 

to  the 

number 

of  lone  pairs  on  the  central  atom. 
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Table  4.29.  Experimental  and  calculated  (atom  charge  model) 


relaxation 

energy  shifts 

(eV)  -  S 

.  b 
series 

• 

Exptl . 

CNDO/2 

EWMO 

Molecule 

AR+AR  * 

.Li  JLi 

ARt 

JLt 

iRL* 

arl 

arl* 

H2S 

0 

0 

0 

0 

0 

(CH3)  sh 

1.31 

1.05 

0.74 

0.69 

0.63 

(ch3)2s 

2.26 

1.81 

1.50 

1.28 

1.17 

CH3SSCH3 

3.15 

2.43 

1.81 

1.68 

1.64 

CF3SSCF3 

3.15 

2.75 

2.10 

2.15 

2.14 

SF4 

0 

0 

0 

0 

0 

S°F2 

-0.56 

-0.76 

-1.02 

-0.20 

-0.33 

(CH3) 2S0 

0.93 

0.76 

-0.71 

0.52 

0.13 

S0C12 

1.75 

1.21 

0.44 

0.78 

0.65 

S02 

-0.90 

-1.57 

-1.91 

-0.27 

-0.58 

(CH30) 2S0 

1.26 

1.26 

0.17 

0.41 

0.48 

SF6 

0 

0 

0 

0 

0 

S02C1F 

-0.11 

0.01 

-0.56 

0.05 

-0.12 

so2ci2 

0.83 

0.77 

0.28 

0.45 

0.28 

S02F2 

-0.76 

-0.98 

-1.22 

-0.38 

-0.58 

SF  Cl 

D 

1.04 

0.55 

0.44 

0.36 

0.40 

(CH30)2S02 

0.71 

1.19 

0.16 

0.21 

0.13 

(a)  Theoretical  values 

for 

k  used  (Table  4.10) 

• 

(b)  The  molecules  have 

been 

i  grouped 

according 

to  the 

number 

of  lone  pairs  on  the  central  atom. 
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comparison  of  similar  molecules  containing  F  and  Cl 
substituents,  or  comparison  of  molecules  with  a  terminal 
sulphur  with  those  containing  the  more  electronegative 
terminal  0  in  the  phosphorus  series  or  comparison  of  the 
-PF^  ligand  with  the  -PCl^  in  the  sulphur  series. 

Other  factors,  however,  appear  to  be  equally  impor¬ 
tant.  Based  upon  the  simple  electronegativity  arguments 

it  might  be  expected  that  -C_H_  and  -CH_  would  behave  sim- 

Z  D  j 

ilarly  and  that  PF^  would  have  substantially  less  re¬ 
laxation  than  PH^/  as  would  molecules  with  a  -CF^  group 
when  compared  to  similar  molecules  containing  a  -CH^ 
group.  This  is  not  the  case,  the  group  exhibits 

more  relaxation  than  the  -CH^  group,  and  in  the  latter  two 
examples  the  reverse  is  observed.  These  can  be  attributed 
to  an  increase  in  the  polarisability  of  the  group  due  to 
the  presence  of  a  larger  number  of  electrons  thereby  con¬ 
tributing  a  larger  total  "electron  flow". 

Replacement  of  a  lone  pair  by  two  F  groups  (e.g.  com¬ 
paring,  respectively,  PF^  SF^,  SC^  with  PF,.,  SF^,  S°2F2^ 
increases  the  number  of  electrons  able  to  "flow"  onto  the 
ionised  centre  and  hence  increases  the  relaxation  in  spite 
of  the  high  electronegativity  of  the  F  groups.  This  is  also 
the  case  when  replacing  the  lone  pair  with  a  terminal  S 
group,  however,  substitution  of  a  terminal  0  group  for  the 
lone  pair  produces  little  net  effect;  presumably  the  electro- 
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negativity  of  0  cancels  out  the  increase  in  relaxation  due 
to  the  presence  of  more  electrons. 

Obviously  relaxation  is  the  sum  of  many  complex  con¬ 
tributions.  The  above  examples  have  been  discussed  in  terms 
of  "number  of  electrons"  and  electronegativity.  Such  factors 
as  "single"  and  "double"  bonds  allowing  more  or  less 
electron  flow  could  also  contribute  to  the  total  amount  of 
relaxation. 

One  final  interesting  point  to  note  is  the  reversal  of 
relaxation  effects  between  -OCH^  and  -CH^  in  the  comparison 
of  the  phosphorus  and  sulphur  series  of  compounds.  The 
EWMO  (Huckel)  calculation  correctly  predicts  this.  The 
calculation  shows  that  in  PCOCH^)^  and  PCCH^)^  the  dif¬ 
ference  in  relaxation  between  them  is  approximately  the 
same  for  both  ionisation  steps  whereas  in  (CH^O^SO  there 
is  an  increase  in  the  relaxation  for  the  second  ionisation 
compared  to  that  for  (CH^^SO. 


■  ; 
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Chapter  5.  Summary  and  Conclusions 


The  KL^L^  Auger  line.  Is  and  2p  core  electron 
lines  were  collected  in  a  series  of  phosphorus  and 
sulphur  compounds  and  their  chemical  shifts  were 
compared.  The  potential  model  was  used  to  attempt  to 
explain  the  shifts. 

The  correlation  between  the  Is  photoelectron  shift 
and  the  2p  photoelectron  shift  was  very  good.  The 
Is  photoelectron  shift  was  larger.  There  was  no  such 
correlation  when  comparing  the  Auger  electron  shift  to 
the  core  electron  shifts.  The  Auger  parameter, ^ 
when  defined  properly,  related  the  Auger  electron  shift 
to  the  core  electron  shift  in  a  convenient  way  and  the 
variations  in  shift  could  be  rationalised  by  con¬ 
sideration  of  the  shifts  in  relaxation  due  to  the  loss 
of  two  electrons.  It  was  found  that  the  relaxation  de¬ 
pends  upon  the  ability  of  a  group  attached  to  the  central 
atom  to  release  electrons.  Many  factors  can  contribute 
to  this  and  it  is  necessary  to  keep  all  of  them  in  mind. 
For  instance,  the  results  showed  that  an  increase  in 
the  polarisability  of  a  group  due  to  the  presence  of  a 
large  number  of  electrons  can  reverse  the  expected  re¬ 
laxation  trends  based  solely  upon  electronegativity. 

The  lack  of  success  of  the  group  shift  concept  when 
related  to  a  large  series  of  molecules  indicated  that 
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the  groups  are  affected  by  the  other  substituents  on  the 
central  atom  and  that  they  cannot  be  ignored. 

The  models  were  tested  with  the  aid  of  two  types  of 
semi-empirical  methods.  The  CNDO/2  calculation  gave 
reasonable  correlation  when  using  fitted  parameters 
though  the  values  obtained  for  k  (the  coulombic  repulsion 
integral)  did  not  match  the  expected  theoretical  trends  too 
well.  It  was  thus  difficult  to  relate  the  relaxation  as¬ 
sociated  with  the  loss  of  the  second  electron  to  that  as¬ 
sociated  with  the  initial  photoionisation.  The  EWMO  (Huckel) 
calculation  was  found  to  partion  charge  differently  depending 
upon  the  number  of  lone  pairs  present  on  the  ionised  atom 
and  thus  the  general  applicability  of  this  method  is 
limited . 

It  is  apparent  that  more  work  is  required  in  form¬ 
ulating  semi-empirical  approaches  to  give  "correct" 
charge  distributions,  especially  if  one  wants  to  relate 
the  relaxation  processes  involved  with  the  ions  of  the 
Auger  electron  to  the  relaxation  processes  involved  with 
the  loss  of  initial  core  electron. 
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APPENDIX  1 


Molecular  Geometries  Used  in  the  Calculations 

Where  experimental  geometries  were  available  only 
the  reference  is  given.  These  are  presented  in  Table  Al.l. 
Where  the  bond  lengths  and  angles  had  to  be  estimated, 
the  assumptions  made  are  also  presented  (Table  A1.2). 
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Table  Al.l.  References  from  which  experimental  geometries 


were  obtained. 

Reference3 

Compound 

Al.l 

ph3. 

PF3,  PF5,  OPF3,  0PC13,  SPC13,  (CH30)3P0, 

sci2, 

sf4,  SiCl4 

A1.2 

sv 

SF5C1,  S0C12,  S02CX2 

A1.3 

H2S' 

(CH3)SH,  (CH3)2S,  (CH3)2SO,  (CH3)SS(CH3) 

sof2, 

S02F2,  S02C1F 

A1.4 

P(CH3 

*  3 

A1.5 

p(cf3 

*  3 

A1.6 

spf3. 

PC13,  (CF3)SS(CF3) 

A1.7 

so2 

aReferences  are: 

Al.l  Landolt-Bornstein  (New  Series)  Group  II,  Vol.  7, 
Springer-Verlag ,  Berlin  (1976) . 

A1.2  P.H.  Laur  in  "Sulphur  in  Organic  and  Inorganic  Chem¬ 
istry",  (ed.  A.  Senning) ,  Vol.  3,  Chap.  24,  Marcel 
Dekker,  New  York  (1972) . 

A1.3  M.D.  Harmony,  V.W.  Laurie,  R.L.  Kuczkowski,  R.H. 

Schwendeman,  D . A .  Ramsay,  F.J.  Lovas,  W.J.  Lafferty 
and  A . G .  Maki,  J.  Phys.  Chem.  Ref.  Data,  8,  619  (1979). 


Cont ' d . 
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Table  Al.l.  Cont'd. 


Al.4  I.M.  Hillier  and  V.R.  Saunders,  Trans.  Far.  Soc.,  66 , 
2401  (1970). 

Al.5  C.J.  Marsden  and  L.S.  Bartlett,  Inorg.  Chem. ,  15 , 

2713  (1976). 

Al.6  "Tables  of  Interatomic  Distances  and  Configuration 
in  Molecules  and  Ions",  (ed.  L.E.  Sutton),  Chem. 

Soc.  Spec.  Publ.  #11  (1958)  . 

A1.7  L.  Asplund,  P.  Kelfve,  B.  Blomster,  H.  Siegbahn, 

K.  Siegbahn,  R.L.  Lozes  and  U.I.  Wahlgren,  Physica 
Spectra,  16,  273  (1977). 
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APPENDIX  2 


Reparameterisation  of  the  EWMO  (Huckel)  Program 

The  program  was  a  modified  iterative  extended  Huckel 
39 

program  (EWMO) .  The  major  features  of  this  program  have 

24  48 

been  adequately  described  '  and  will  not  be_dealt  with 

here  with  the  exception  of  the  parameterisation . 

The  core-energy  parameters  and  the  values  for  the 

electron  repulsion  integrals  were  taken  from  Sichel  and 
49 

Whitehead.  These  atom  parameters  are  based  upon  valence 

state  energies  for  the  atom  given  by  Hinze  and  Jaffe.^ 

Figure  A2 . 1  shows  a  plot  of  the  core-integrals,  U  and 

s  s 

U  ,  against  the  second-row  elements.  The  manner  in  which 
PP 

the  argon  values  were  obtained  has  already  been  described 
in  chapter  3 . 

The  graph  shows  a  smooth  trend  of  values  with  the 

exception  of  those  pertaining  to  phosphorus.  A  similar 

57 

plot  for  the  values  obtained  by  Oleari  et  al.  is  shown  in 
Figure  A2 . 2 .  In  this  case  a  smooth  curve  through  all  the 
points  is  obtained. 

As  mentioned  previously  in  chapter  3,  the  electron 

repulsion  integrals,  g  ,  in  the  CNDO  formalism  become 

proportional  to  the  Slater  exponent,  Cg,  and  so  a  plot  of 

a  versus  l  for  the  second  row  should  be  linear.  Sichel 
yss 

and  Whitehead, 49  in  their  parameterisation,  obtained  values 
for  an  average  electron  repulsion  integral,  g^,  for  the 
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Figure  A2 . 1  Core  integrals  (-U  „  and  -U  )  taken  from 
— 2 -  J  s  s  pp 

49 

Sichel  and  Whitehead  versus  the  second  row  elements. 


-U/eV 


Figure  A2 . 2  Core  integrals  (-Uf,r  and  -U^)  taken  from 

- — — — — ■ — -  O  O  U  M 


57 


versus  the  second  row  elements 


Oleari  et  al. 
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* 

one-centre  term  and  g^  for  the  two-centre  term.  These 
values  are  defined  as 


AA 

1 

28 

(g  + 
^ss 

12g 

^sp 

+  3g 

PP 

+  12g  ,) 

PP 

* 

AA 

1 

16 

(g  + 

^  ss 

6g  + 

^sp 

3g  + 

PP 

6g  i ) 
^pp' 

where  are  the  appropriate  electron-repulsion  integrals 

between  the  s  and  p  electrons.  It  is  not  unreasonable  to 

assume,  to  the  level  of  the  semi-empirical  approximation 

used  here,  that  the  average  electron  repulsion  terms  are 

also  proportional  to  the  Slater  exponent.  A  plot  of  g 

versus  Cg  is  shown  in  Figure  A2 . 3  for  both  sets  of  para- 

57 

meters.  The  parameter  values  given  by  Oleari  et  al.  are 

approximately  linear  for  all  elements  in  the  row  whereas 

49 

in  the  values  given  by  Sichel  and  Whitehead  the  value 
assigned  to  phosphorus  deviates  substantially  from  the 
linearity  of  the  parameters  assigned  to  the  remaining 
elements.  This  deviation  of  the  phosphorus  parameter  in 
the  latter  case  seems  to  reflect  the  stability  of  the  half- 
filled  shell  in  atomic  phosphorus.  However,  when  dealing 
with  changes  in  the  molecular  environment  the  concept  of 
the  stability  of  the  half-filled  shell  is  no  longer  relevant 
because  all  systems  are  based  on  filled  valence  shells. 

Thus  for  the  series  from  Si  through  to  Cl  a  smooth,  gradual 
increase  in  parameter  values  would  be  expected.  In  this  con¬ 
text  the  trends  exhibited  by  the  parameterisation  of  Oleari 
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Slater  Exponent  ( C  3 s ) 


Figure  A2 . 3  Average  one-centre  repulsion  integral  (g  ) 
versus  the  Slater  orbital  exponent  (C3g)  of  the  second  row 
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et  al.^7  seems  more  reasonable. 

Rather  than  attempt  to  reparameterise  the  entire  program 
using  the  values  given  by  Oleari  et  al.,57  we  thought  that 
it  would  be  more  efficient  and  equally  justified  to  inter¬ 
polate  new  values  for  the  phosphorus  parameters,  in  effect 

obtaining  new  parameters  for  phosphorus  which  would  be  related 

49 

to  the  other  atomic  parameters  given  by  Sichel  and  Whitehead. 

57 

This  set  in  general  would  parallel  those  of  Oleari  et  al. 

The  one-centre  electron  repulsion  integral  was  estimated  from 
a  least-squares  fit  of  atomic  parameters  for  the  second 
row  (with  the  exception  of  P)  as  shown  in  Figure  A2 . 3 . 

A  similar  procedure  was  used  to  estimate  the  two-centre 
electron  repulsion  integral.  A  correlation  of  the  core 
integrals  from  both  works  is  shown  in  Figure  A2.4.  The 
phosphorus  parameter  value  given  by  Oleari  et  al.  was  then 
converted  to  an  interpolated  "Sichel  and  Whitehead"  core 
integral  value  for  phosphorus  using  the  least-squares 
correlation  function  defined  by  the  data  in  Figure  A2 . 4 . 

The  final  values  for  the  new  phosphorus  parameters  are 
presented  in  Table  A2 . 1 .  Figure  A2 . 5  shows  the  new  relation¬ 
ship  of  the  core  integrals  for  the  second  row  elements. 


•V  . 
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Figure  A2 . 4  Core  integrals  (-U  and  -U  )  of  Oleari  et 

49 

versus  those  of  Sichel  and  Whitehead. 
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Table  A2.1 


New  Parameters  for  P  (in  eV) 


-U 

ss 


P  50.648  42.949 


8.079 


8.124 


.  $  5  ^9IH I  ’  H 
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Figure  A2 . 5  Core  integrals  (-Ugs  and  -U  )  taken  from  Sichel 
and  Whitehead49  with  the  new  P  parameterisation  versus  the 


second  row  elements. 
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APPENDIX  3 


Corrected  values  for  the  calculated  Auger  shifts  as 
presented  in  Table  1,  ref.  8. 

Auger 

2  step  3  step 

egn.  5  eqn.  6 


PH3 

0 

0 

PC13 

-1.78 

-1.32 

PF3 

5.08 

5.53 

SPCI3 

0.43 

0.60 

spf3 

5.90 

6.33 

OPCI3 

0.72 

1.13 

OPE  3 

7.54 

7.96 

PFC 

D 

8.64 

9.28 
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